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ABSTRACT 


The  ability  to  produce  an  accurate  estimate  of  the  expected  amount  of  damage  that 
will  occur  to  a  targeted  region  from  the  fragments  generated  by  intercepting  a  Theater/ 
Tactical  Ballistic  Missile  (TBM)  in  flight  is  necessary  for  determining  whether  the  TBM 
should  or  should  not  be  intercepted.  For  a  sea-based  Theater  Missile  Defense  (TMD) 
system  that  is  capable  of  predicting  the  remaining  portion  of  the  TBM's  trajectory  along 
with  its  impact  point  from  its  current  and  past  flight  parameters,  possessing  this  ability 
may  prevent  unnecessary  intercepts,  and  thus  unnecessary  damage.  This  thesis  proposes  a 
simplistic  simulation  model  that  produces  expected  damage  estimates  for  various  intercept 
ranges.  Damage  is  evaluated  in  three  areas;  number  of  fragments  to  impact  the  target 
region,  total  mass  to  impact  the  target  region,  and  total  kinetic  energy  to  impact  the  target 
region.  The  values  generated  by  this  model  can  aid  in  the  decision  of  whether  to  intercept 
the  TBM  or  not  by  comparing  them  against  the  expected  damage  values  caused  by  the 
unintercepted  TBM  which  are  calculated  by  the  sea-based  TMD  system. 
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EXECUTIVE  SUMMARY 


PROBLEM:  A  model  for  estimating  the  expected  collateral  damage  to  a  target  region 
resulting  from  the  intercept  of  a  Theater/Tactical  Ballistic  Missile  (TBM)  is  necessary  for 
the  decision  maker  to  determine  if  intercepting  the  TBM  is  the  proper  course  of  action 
and  is  the  focus  of  this  thesis. 

With  the  increased  threat  of  theater/tactical  ballistic  missile  attacks  against 
defenseless  regions,  and  the  propensity  of  some  nations  to  use  ballistic  missiles  as  a 
bullying  tactic,  theater  missile  defense  (TMD)  has  become  a  high-priority  issue  for  the 
Department  of  Defense.  Presently,  IBM's  pose  a  great  threat  to  U.S.  military  operations 
because  almost  every  conceivable  future  conflict  involving  U.S.  military  deployments 
would  involve  the  risk  of  TBM  attacks  against  U.S.  forces.  The  only  existing  deployable 
TMD  system  which  intercepts  a  TBM  in  flight.  Patriot,  an  Army  air  defense  system 
adapted  for  the  TMD  role,  has  suffered  severe  criticism  for  its  performance  in  Desert 
Storm.  At  first.  Patriot  was  praised  as  an  extremely  effective  interceptor  of  SCUD 
missiles  launched  from  Iraq  at  Israel  and  Saudi  Arabia.  Unfortunately,  post-war  analysis 
showed  that  the  Patriot  missile  was  much  less  effective  than  was  initially  claimed.  Patriot 
is  now  being  modified  to  enhance  its  capabilities  against  TBM's.  Another  land-based  TMD 
system  under  development  is  the  new  Theater  High  Altitude  Area  Defense  system 
(THAAD).  The  integration  of  THAAD  and  an  improved  Patriot  promise  to  provide 
outstanding  defense  against  fUture  ballistic  missiles. 

In  addition  to  these  systems,  the  Ballistic  Missile  Defense  Organization  (BMDO) 
has  directed  the  Navy  to  adopt  this  two-tiered  approach  of  exo-atmospheric  and 
endo-atmospheric  defense  and  develop  its  own  TMD  systems  to  be  integrated  with  the 
Army  TMD  systems  in  the  future.  The  weapon  chosen  for  the  lower-tier,  or  endo- 
atmospheric,  role  is  SM-2  Block  IV  A  surface-to-air  missile.  One  concern  with  the  future 
testing  of  this  missile  is  the  unknown  magnitude  of  the  collateral  damage  that  results  from 
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a  TBM  intercept.  Therefore,  a  model  is  proposed  for  estimating  the  expected  damage  to 
a  target  region  and  its  surrounding  areas. 

SOLUTION  METHODOLOGY :  This  thesis  addresses  the  question  of  when  to 
intercept  a  TBM  and  when  not  to  intercept  it.  The  model  proposed  is  a  simplistic 
analytical  attempt  at  estimating  the  distribution  of  the  expected  damage  that  occurs  to  a 
target  region  and  its  surrounding  areas  from  intercepting  a  TBM  at  predetermined  points 
on  its  flight  path.  The  estimate  can  then  be  compared  against  the  amount  of  damage  that 
would  have  occurred  if  the  TBM  were  allowed  to  continue  without  intercept,  assuming 
the  TMD  system  in  question  has  this  capability. 

Utilizing  MODSIM,  a  simulation,  which  replicated  the  behavior  of  the  fragments 
generated  after  an  intercept,  was  developed.  Because  of  the  lack  of  available  data 
concerning  the  true  fragmentation  of  a  TBM  upon  intercept,  many  simplifying  assumptions 
were  made.  In  order  to  account  for  the  uncertainty  in  the  TBM's  impact  point,  if  not 
intercepted,  two  parameters  of  the  TBM's  initial  conditions  were  allowed  to  vary  in  each 
run:  vertical  launch  angle  and  off-target  launch  angle.  Replications  of  the  simulation  were 
conducted  for  various  intercept  points  and  pairs  of  launch  angles  so  as  to  generate  data  on 
the  number  of  hits,  total  mass  and  total  kinetic  energy  impacting  the  target  region.  The 
distributions  of  the  number  of  fragment  hits,  the  total  mass  impacting  the  target  region  and 
the  total  kinetic  energy  impacting  the  target  region  are  examined  to  determine  expected 
values. 

SOLUTION:  The  modeling  and  simulation  results  presented  are  speculative  and 
exploratory.  They  illustrate  the  effects  of  various  physical  influences  upon  missile 
collateral  damage  to  a  target  in  a  preliminary  way.  The  conclusion  drawn  from  this  thesis 
is  that  knowing  the  expected  damage  to  occur  to  a  target  region  if  a  TBM  is  intercepted  is 
helpful  in  making  a  decision  whether  to  intercept  it  or  not.  This  model  is  a  simple 
analytical  attempt  to  produce  a  tentative  estimate.  At  this  time,  very  little  is  known 
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concerning  the  dispersion  of  fragments  generated  by  an  intercept.  Because  of  this  lack  of 
information,  a  crude  estimate  is  the  only  estimate  available.  When  the  assumptions  made 
in  this  thesis  are  determined  with  certainty,  a  more  detailed  version  of  the  model, 
replicating  the  exact  flight  paths  and  impact  points,  can  be  developed. 
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1.  INTRODUCTION 


With  the  increased  threat  of  theater/tactical  ballistic  missile  (TBM)  attacks  against 
defenseless  regions,  and  the  propensity  of  some  nations  to  use  ballistic  missiles  as  a 
bullying  tactic,  theater  missile  defense  (TMD)  has  become  a  high-priority  issue  for  the 
Department  of  Defense.  Presently,  TBM's  pose  a  great  threat  to  U.S.  military  operations 
because  no  truly  effective  ballistic  missile  defense  systems  exist.  Moreover,  almost  every 
possible  future  conflict  involving  U.S.  military  deployments  would  also  involve  the  risk  of 
TBM  attacks  against  U.S.  forces.  The  most  familiar  TBM,  encountered  extensively 
during  the  Gulf  War,  is  the  SCUD,  Variant  B.  Many  countries  have  these  missiles  in  their 
arsenals  including  such  unstable  nations  as:  North  Korea,  Iraq,  Vietnam,  Yemen,  Libya, 
Afghanistan,  and  Syria  [Ref  1]. 

The  only  existing  TMD  system  which  intercepts  a  TBM  in  flight,  the  Patriot,  an 
Army  air  defense  system  adapted  for  the  TMD  role,  has  suffered  severe  criticism  for  its 
performance  in  Desert  Storm.  At  first.  Patriot  was  praised  as  an  extremely  effective 
interceptor  of  SCUD  missiles  launched  from  Iraq  at  Israel  and  Saudi  Arabia. 

Unfortunately,  post-war  analysis  showed  that  the  Patriot  missile  was  much  less  effective 
than  was  initially  claimed.  In  fact,  the  majority  of  SCUDs  thought  to  be  intercepted  had 
spontaneously  broken  into  segments  upon  re-entry  as  a  result  of  faulty  manufacturing. 
Patriot  is  now  being  modified  to  enhance  its  capabilities  against  TBM's.  Another 
land-based  TMD  system  under  development  is  the  new  Theater  High  Altitude  Area 
Defense  system  (THAAD).  The  integration  of  THAAD  and  an  improved  Patriot  promise 
to  provide  outstanding  defense  against  future  ballistic  missiles.  However,  THAAD 
development  is  being  decelerated  because  of  budget  cuts.  The  drawback  of  both  Patriot 
and  THAAD,  when  the  latter  is  completed,  is  the  enormous  lead  time  required  to  deploy 
and  set  up  these  systems.  Future  TBM  attacks  on  areas  of  U.S.  interest  may  not  be 
known  weeks  in  advance;  time  to  deploy  these  defense  systems  may  not  be  available.  Also, 
permission  must  be  granted  from  foreign  leaders  to  deploy  land-based  TMD  systems.  The 
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leaders  may  be  hesitant  even  though  these  systems  are  there  to  save  their  territory  because 
deployment  of  such  systems  on  their  land  may  escalate  hostilities.  In  addition  to  the 
deployment  lead  time  is  the  enormous  expense  of  airlifting  the  systems  to  the  threat  areas. 
In  summary,  such  drawbacks  provide  strong  motivation  for  the  development  of  a 
sea-based  ballistic  missile  defense  system. 

As  mentioned  previously,  plans  and  initial  steps  necessary  for  TBM  attacks  may 
not  be  known  long  enough  (e  g.,  weeks)  in  advance  of  attack  to  allow  for  deployment  of 
Patriot  batteries.  A  timely  response  option  is  the  use  of  a  mobile  sea-based  ballistic  missile 
defense  system.  The  Navy  is  in  the  process  of  developing  a  version  of  this  system  by 
upgrading  the  Aegis  Weapon  System  and  the  SM-2  Block  IV  standard  missile  in  order  to 
effectively  detect,  track,  engage,  and  intercept  TBM's.  The  new  standard  missile  is  to  be 
an  improvement  of  the  existing  vertical-launch  Terrier/standard  surface-to-air  missile.  The 
latter  requires  modifications  of  the  warhead,  fuze,  and  seeker  to  improve  its  capability  to 
defend  against  TBM's.  The  new  version  would  build  upon  already-developed  engineering 
and  design  efforts  used  on  Patriot,  but  there  are  drawbacks  to  using  this  technology  [Ref 
1].  One  of  these  is  the  explosive  warhead.  This  warhead  may  cause  total  destruction  with 
a  near-miss  intercept,  or  it  may  not.  If  the  near-miss  intercept  does  not  totally  destroy  the 
TBM,  debris  from  the  TBM  could  impact  the  intended  target  or  surrounding  areas. 

Collateral  damage  is  the  term  for  non-target  damage  caused  by  the  debris.  For 
example,  if  a  TBM  is  targeted  to  impact  the  center  of  a  coastal  city  and  U.S.  forces  are 
conducting  an  amphibious  assault,  intercept  of  the  TBM  might  stop  it  from  directly 
impacting  its  target,  but  the  fragments  generated  by  the  intercept  might  fall  on  the  U.S. 
forces  and  interfere  with  the  assault  operation.  When  contemplating  an  intercept,  the 
pattern  of  the  debris  generated  from  the  intercept  must  be  estimated  in  order  to  determine 
the  amount  of  damage  that  will  occur  if  the  ballistic  missile  is  intercepted  and  fragments 
continue  on  course.  According  to  Israeli  reports  after  Desert  Storm,  the  1 1  SCUDS 
intercepted  by  Patriot  batteries  caused  more  damage  than  did  the  13  SCUDS  that  directly 
hit  targets  within  Israel.  This  raises  an  interesting  question:  with  prior  knowledge  as  to 
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the  amount  of  damage,  both  primary  and  collateral,  that  will  occur  from  intercepting  a 
TBM,  is  intercepting  the  TBM  always  the  best  strategy?  If  the  TBM  is  allowed  to 
continue  on  course  without  intercept,  it  is  possible  that  the  amount  of  damage  resulting 
from  the  single  impact  point  may  be  less  than  that  of  a  shower  of  fragments  impacting  the 
target  area.  Of  course,  there  is  also  a  chance  that  the  TBM  will  miss  the  target  area 
completely.  Because  the  dispersion  of  fragments  resulting  from  an  intercept  and  the  size 
of  each  fragment  are  unknown,  an  accurate  estimate  of  the  debris  pattern  is  almost 
impossible.  This  dispersion  uncertainty  may  be  caused  by  numerous  aspects  including 
angle  of  impact  by  interceptor,  angle  of  descent  of  ballistic  missile,  impact  point  on 
ballistic  missile,  distance  between  interceptor  and  ballistic  missile  at  intercept. 

Purpose  of  this  Thesis:  Where,  and  If,  to  Intercept? 

In  order  to  address  the  question  as  to  when  to  intercept  a  ballistic  missile  and  when 
not  to  intercept  it,  a  simulation  model  is  developed  to  estimate  the  distribution  of  the 
damage  that  occurs  from  intercepting  a  TBM  at  predetermined  points  on  its  trajectory. 

This  damage  is  then  compared  to  the  amount  of  damage  that  would  have  occurred  if  the 
TBM  were  allowed  to  continue  without  intercept.  Development  of  the  model  requires  a 
number  of  simplifying  assumptions  to  be  made.  The  first,  and  possibly  most  credible,  set  of 
assumptions  deals  with  the  characteristics  of  the  TBM,  including  launch  acceleration, 
motor  burn  time,  missile  length,  missile  diameter  and  missile  mass.  The  second,  more 
speculative,  set  of  assumptions  concerns  the  generation  and  motion  of  fragments, 
including  number  of  fragments,  mass  of  each  fragment,  cross-sectional  area  of  each 
fragment,  and  the  velocity  of  each  fragment  after  intercept.  These  characteristics  are 
generated  from  various  probability  distributions  judged  to  be  appropriate.  The  final 
assumption  is  target  region  designation.  For  this  study  the  target  region  is  a  coastal  region 
that  is  2  km  X  2  km.  This  region  is  then  broken  up  into  16  blocks  that  are  0.5  km  X  0.5 
km  each.  Actual  targets  that  are  subregions  of  this  region  can  be  studied  for  damage. 
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The  preceding  assumptions  remain  the  same  for  each  run  of  the  model.  However, 
two  parameters  are  allowed  to  vary  in  each  run:  vertical  launch  angle  and  off-target 
launch  angle.  For  each  pair  of  angles,  while  holding  all  other  parameters  constant,  a 
different  IBM  flight  path  is  generated.  From  these  flight  paths  an  impact  point  could  be 
calculated.  Although,  having  precise  knowledge  of  these  two  angles  is  not  realistic, 
determining  the  impact  point  of  a  IBM  from  its  flight  path,  conditional  on  the  angles,  is 
realistic.  For  simulation  purposes  the  two  angles  are  taken  to  be  given  and  known  in  order 
to  calculate  the  impact  point.  Realistic  expected  impact  damage  can  be  estimated  by 
averaging  over  any  angle  distribution  of  interest. 

The  model  is  intended  to  represent  a  worst-case  situation:  the  intercept  results  in  a 
near-miss  warhead  explosion  that  creates  missile  fragments  which  continue  traveling 
towards  the  target  area  with  minimal  loss  of  mass  and  minimal  change  in  velocity.  The 
TBM's  flight  path  is  modeled  by  using  suitable  equations  of  motion  in  three  dimensions. 

At  pre-determined  points  on  the  flight  path,  the  hypothetical  intercept  points,  the  missile 
breaks  up  into  a  random  number  of  fragments  resulting  from  the  near-miss  intercept.  Each 
fragment  has  the  velocity  of  the  TBM  at  intercept  plus  a  random  "kick"  in  three 
dimensions  to  simulate  breaking  apart.  Each  of  the  fragments  continues  towards  the 
target  region  on  a  path  determined  by  its  mass  and  cross-sectional  area.  This  model  treats 
the  cross-sectional  areas  of  the  fragments  as  functions  of  their  respective  masses,  i.e.,  the 
larger  the  mass  the  larger  the  cross-sectional  area. 

The  objective  of  this  modeling  project  is  to  estimate  the  expected  damage  that 
occurs  when  the  TBM  is  intercepted  at  various  distances  from  the  target  area,  given  the 
vertical  and  off-angle  launch  angles.  Certain  intercept  points  generate  fragments  which 
impact  the  target  area.  The  amount  of  damage  that  occurs  from  a  fragment  impact  is 
estimated  using  a  damage  function  that  takes  into  account  fragment  mass  and  velocity. 

The  damage  generated  from  different  intercept  points  is  then  compared  to  the  damage 
generated  by  letting  the  missile  through  the  defenses  with  hope  of  a  miss.  Replications  of 
this  simulation  are  conducted  for  various  intercept  points  and  pairs  of  launch  angles  so  as 
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to  generate  data  on  the  numbers  of  hits,  total  mass  and  total  kinetic  energy  impacting  the 
target  area.  The  distributions  of  the  number  of  fragment  hits,  the  total  mass  impacting  the 
target  area  and  the  total  kinetic  energy  impacting  the  target  area  will  be  examined  to 
determine  expected  values  and  summaries  of  the  effect  of  intercepting  missiles  at  various 
locations  along  their  trajectories.  The  ballistic  missile  modeled  in  these  simulations  are 
based  on  the  characteristics  of  the  SCUD  variants  used  in  Desert  Storm,  i.e.,  they  will 
have  a  fairly  large  circular  error  probable  (CEP)  which  is  represented  by  the  varying 
angles 

The  modeling  and  simulation  results  presented  are  speculative  and  exploratory. 
They  illustrate  the  effects  of  various  physical  influences  upon  missile  collateral  damage  to 
a  target  in  a  preliminary  way.  The  model  equations  and  software,  which  are  provided 
herein,  can  be  modified  or  extended  to  test  sensitivity  to  specific  parameter  values.  With 
slightly  more  effort  the  actual  physical  relationships  that  have  been  used  can,  if  considered 
necessary,  be  modified  in  more  reasonable  and  defensible  ways. 
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II.  NATURE  OF  THE  PROBLEM 


A.  PROBLEM  DESCRIPTION 

As  demonstrated  in  the  Gulf  War,  Iraq  used  it's  A1  Hussein  TBMs,  i.e.  modified 
SCUD  B's,  against  civilian  targets,  most  of  which  were  populated  cities.  The  circular  error 
probable  (CEP)  of  this  particular  TBM  is  fairly  large,  being  about  1000  m,  allowing  a 
chance  of  missing  the  target  altogether  and  landing  in  some  insignificant  area,  e.g.,  the 
sea.  Actually,  this  happened  quite  often  during  the  Gulf  War.  Of  the  88  TBMs  launched 
by  Iraq  towards  Israel  and  Saudi  Arabia,  35  missed  their  targets  completely,  landing  either 
in  the  sea  or  in  the  desert.  Future  sea-based  TMD  systems  will  be  able  to  estimate,  with  a 
level  of  certainty,  the  impact  point  of  a  TBM  by  calculating  the  remaining  portion  of  its 
trajectory  according  to  the  TBM's  current  and  past  parameters.  However,  these  systems 
do  not  rely  solely  on  the  estimate  in  determining  what  action  to  take,  since  the  estimate  is 
not  exact.  If  the  impact  point  is  judged  likely  to  be  close  to  the  target  area,  one  of  two 
courses  of  action  must  be  taken:  intercept  the  missile  and  risk  collateral  damage  by  the 
missile  fragments,  or  do  not  intercept  the  missile  and  hope  nothing  significant  is  located 
near  the  impact  point.  Therefore,  the  problem  to  be  addressed  in  this  section  is  to 
determine  whether  a  sea-based  TMD  system,  assigned  to  protect  a  targeted  region,  e.g.,  a 
coastal  city,  should  intercept  a  TBM  in  flight,  creating  missile  fragments  that  may  still 
impact  the  target  area,  or  should  refrain  from  interception,  allowing  the  TBM  to  continue 
towards  the  intended  target  area,  having  assessed  that  it  will  miss. 

To  address  this  problem,  a  computer  simulation  utilizing  MODSIM,  a  simulation 
software  system,  is  designed  to  determine  the  amount  of  damage  that  occurs  from 
intercepting  a  TBM  in  flight.  In  this  simulation  a  TBM  modeled  after  the  A1  Hussein,  but 
with  generic  numbers  substituted  for  some  of  its  characteristics  to  preserve  the 
unclassified  status  of  this  project,  is  launched  toward  a  predefined  target  area.  Parameters 
of  the  TBM  are  first  chosen  to  allow  an  impact  point  at  the  exact  center  of  the  target  area. 
Points  along  the  TBM's  flight  path  are  selected  as  intercept  points,  at  which  the  missile 
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separates  into  a  random  number  of  fragments  that  continue  towards  the  target.  The  final 
position,  or  impact  point,  of  each  of  the  fragments  generated  is  examined  to  determine  if 
the  fragment  landed  inside  the  targeted  area,  and  where  therein.  Summing  the  results  for 
each  fragment  totals  the  amount  of  damage  to  occur  for  a  given  intercept  point. 
Conducting  numerous  replications  of  the  model  for  each  intercept  point  yields  the  average 
damage  to  occur  to  the  target  area.  These  results  are  further  analyzed  using  statistical 
techniques  to  determine  the  expected  damage  to  occur  for  various  intercept  points. 

Since  this  model  is  deterministic,  the  unintercepted  TBM,  under  the  assumption  of 
a  flat,  non-rotating  earth,  would  have  the  same  impact  point,  i.e.,  the  center  of  the  target 
area,  for  each  simulation  run  unless  one  or  more  parameters  are  allowed  to  vary.  By 
changing  parameters  on  each  simulation  run,  the  TBM's  impact  point,  if  allowed  to 
continue  towards  the  target  without  intercept,  also  changes.  The  parameters  chosen  to 
vary  in  this  simulation  are  the  vertical  and  off-target  launch  angles  as  depicted  in  Figure  1 . 
By  introducing  a  small  degree  of  error  in  either  angle,  the  TBM  impacts  vary  from  the  aim 
point.  As  mentioned  previously,  a  sea-based  TMD  system  will  be  able  to  estimate,  fairly 
accurately,  this  final  hypothetical  impact  point  of  the  TBM  by  constructing  the  remaining 
trajectory  using  its  current  and  past  parameters;  course,  speed,  angle  of  descent,  altitude, 
and  others.  Hence,  a  decision  can  be  made  either  to  intercept  the  TBM,  and  where  along 
its  path,  or  to  allow  it  to  continue  on  course  by  comparing  the  damage  caused  by  a  TBM 
impact  at  the  point  estimated  by  the  sea-based  TMD  system  against  the  calculated  damage 
caused  by  fragments  impacting  the  target  area. 
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Figure  1 .  3 -Dimensional  Plot  of  the  TBM  Launch  Angles 

t  =  off-target  launch  angle 
V  =  vertical  launch  angle 


B.  PROBLEM  ASSUMPTIONS 

To  preserve  the  unclassified  status  of  this  thesis  and  to  allow  the  expressions 
utilized  in  developing  the  TBM  trajectories  to  be  closed  form,  the  following  assumptions 
are  made: 

-  The  sea-based  TMD  system  is  stationed  somewhere  between  the  launch  point 
and  target  where  alt  potential  intercept  points  utilizing  an  SM-2  Block  IV  A  Standard 
missile  are  feasible.  The  position  of  the  ship  will  not  coincide  with  the  debris  pattern. 


9 


-  Initial  detection  of  the  TBM  is  via  space-based  assets  which  provide  a  cue  to  the 
sea-based  TMD  weapons  system.  The  sea-based  TMD  system  picks  up  the  TBM  track  as 
soon  as  it  enters  the  detection  envelope,  giving  it  enough  time  to  calculate  an  impact 
point.  Error  in  that  estimate  will  realistically  occur,  but  is  not  considered  in  what  follows. 

-  The  attack  on  the  target  consists  of  only  one  TBM. 

-  The  TBM  experiences  constant  acceleration  and  straight-line  motion  during 
boost  phase  with  air  resistance  having  no  effect. 

-  The  TBM  follows  a  normal  ballistic  path,  executing  no  maneuvers  to  avoid  an 
intercepting  missile. 

-  The  TBM  remains  intact  throughout  the  flight  and  the  mass  remains  essentially 
constant  during  boost  phase  even  though  fuel  is  being  expended.  At  the  completion  of 
boost  phase  the  TBM's  mass  is  decreased  by  one  seventh  of  its  original  mass  to 
compensate  for  the  spent  fuel. 

-  The  intercept  simulates  a  near-miss  explosion  of  the  SM-2's  warhead  which 
separates  the  TBM  into  several  fragments,  the  number  and  sizes  are  simulated  from  a 
combination  of  plausible  probability  distributions;  these  latter  can  be  modified  as  desired 
by  an  analyst. 

-  There  is  minimal  loss  of  mass  and  minimal  change  in  velocity  resulting  from  the 
intercept. 

C.  MODEL  OBJECTIVE 

If  it  is  estimated  that  the  TBM  will  impact  in  the  vicinity  of  the  target,  intercepting 
it  towards  the  end  of  its  flight  path  at  low  altitude  may  generate  fragments  that  still  impact 
the  target  region  and  cause  collateral  damage.  Intercepting  it  earlier  in  its  flight  path 
obviously  presents  a  smaller  chance  of  collateral  damage.  With  this  knowledge,  the  ideal 
intercept  is  the  farthest  distance  possible  from  the  target  where  intercept  is  feasible. 
However,  at  this  distance  the  impact  point  of  the  TBM  estimated  by  the  sea-based  TMD 
system  is  fairly  inaccurate,  and  an  unnecessary  intercept  may  occur,  i.e.,  the  TBM  would 
have  fallen  short  of  its  target  and  landed  in  the  sea,  requiring  no  intercept.  Before 
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planning  an  engagement,  knowing  the  damage  that  results  from  different  intercept  points 
allows  the  sea-based  TMD  system  to  decide  whether  an  intercept  is  necessary.  Therefore, 
the  objective  of  this  model  is,  for  particular  intercept  points,  to  estimate  the  expected 
values  for  the  number  of ji'agment  hits,  the  amount  of  mass  to  impact  the  subregions  of 
the  general  target  region,  and  the  total  kinetic  energy  to  impact  subregions  of  the  general 
target  region  in  order  to  decide  whether  the  damage  sustained  from  the  fragments  will  be 
greater  than  that  of  the  TBM  itself  The  above  measures  of  TBM  effect  are  actually 
tabulated  for  subregions  of  the  general  target  region.  Note  that  the  simulation  allows 
evaluation  of  measures  of  target  (sub)region  effects  other  than  expected  values;  it  is,  for 
example,  possible  to  estimate  the  probability  that  the  number  of  fragments  that  impact  in  a 
particular  subregion  exceeds  a  prescribed  level,  or  the  probability  that  the  total  kinetic 
energy  delivered  to  a  particular  subregion  exceeds  a  given  level.  This  capability  of  the 
model  is  appropriate  when  hardened  military  targets  are  in  the  area,  and  their  vulnerability 
to  damage  is  low  unless  a  specific  threshold  is  exceeded. 
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III.  MODEL 


The  purpose  of  this  thesis  is  to  develop  a  simulation  model  that  produces  a 
tentative  estimate  of  the  expected  damage  that  occurs  within  a  designated  target  region 
caused  by  the  fragments  generated  from  intercepting  a  TBM.  To  accomplish  this  task  the 
model  may  be  considered  as  having  two  parts.  The  first  part  represents  launching  a  TBM 
towards  a  designated  target  region  located  some  known  distance  from  a  selected  launch 
point.  The  TBM  is  tracked  from  launch  through  its  trajectory  until  it  reaches  the  intercept 
point  chosen  by  the  user.  The  effect  of  tracking  the  TBM  is  simulated  by  formulating  and 
solving  three-dimensional  equations  of  motion,  utilizing  the  TBM's  mass  and 
cross-sectional  area  and  physical  laws. 

The  second  part  of  the  model  begins  with  the  TBM  intercept,  which  is  treated  as  a 
near-miss  intercept  that  only  breaks  the  TBM  into  pieces  rather  than  destroying  it.  The 
reason  behind  this  is  to  have  the  model  represent  a  worst-case  situation.  If  minimal  mass 
and  velocity  are  lost  as  a  result  of  an  intercept  there  stands  a  greater  possibility  of  the 
fragments  reaching  the  target  region.  At  intercept  the  TBM  breaks  up  into  a  random 
number  of  fragments.  Each  fragment  generated  continues  towards  the  target  region  on  its 
own  flight  path.  These  fragment  flight  paths  are  calculated  using  the  same 
three-dimensional  equations  of  motion  as  describe  the  TBM's  flight  path,  but  this  time 
depending  upon  each  fragment's  mass  and  cross-sectional  area.  Of  concern  in  this  model 
is  the  extent  of  damage  to  the  target  region  caused  by  these  fragments.  Therefore,  the 
model  determines  which  fragments  impact  the  designated  target  region  by  comparing  the 
final  positions  of  each  fragment  to  the  target  location.  If  it  is  determined  that  a  fragment 
has  impacted  the  target  region,  its  mass  and  kinetic  energy  are  recorded  and  subsequently 
totaled  with  the  masses  and  kinetic  energies  of  other  impacting  fragments.  However,  if  it 
is  determined  that  a  fragment  does  not  impact  the  target  region,  it  is  considered  a  miss, 
and  the  mass  and  kinetic  energy  are  recorded  for  that  region  surrounding  the  target  where 
the  miss  occurred. 
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Each  mn  of  the  present  model  is  analyzed  by  statistical  techniques  to  estimate  the 
expected  number  of  fragments  to  impact  various  subregions  of  the  target,  and  the 
expected  total  mass  and  expected  total  kinetic  energy  to  impact  these  subregions. 
Conducting  2000  replications  of  the  model  for  each  intercept  point  enabled  Monte  Carlo 
analysis  to  be  implemented  to  produce  a  fairly  accurate  estimate  of  these  expected  values. 
This  many  replications  were  required  to  keep  the  standard  error  of  the  estimate  at  a 
desired  low  level.  The  following  section  describes  the  model  in  depth. 

A.  SIMULATION  DESCRIPTION 

1.  Initial  Conditions 

As  mentioned  previously,  the  TBM  utilized  for  this  thesis  is  being  modeled  after 
Iraq's  version  of  a  modified  SCUD,  the  A1  Hussein.  For  classification  purposes,  only  four 
actual  characteristics  of  the  A1  Hussein  are  used  in  this  simulation.  The  remaining 
characteristics  required  for  the  simulation  have  been  calculated  based  on  assumptions 
made  concerning  the  distance  to  the  target  region.  Another  assumption  concerning  the 
flight  of  the  TBM  was  that  the  calculations  were  based  on  a  flat,  non-rotating  earth.  For 
this  thesis,  the  only  portion  of  the  TBM's  trajectory  of  interest  is  immediately  pre-intercept 
to  impact.  Therefore,  the  manner  in  which  the  TBM's  flight  path  is  calculated  is  artificially 
simplistic  and  thus  is  performed  as  a  means  of  obtaining  the  trajectory  portion  of  interest. 
The  characteristics  of  the  TBM  modeled  in  this  thesis,  along  with  the  calculations 
necessary  to  compute  the  simplistic  TBM  flight  path  are  contained  in  Appendix  A. 

2.  Intercept  and  Fragment  Generation 

In  order  to  intercept  a  TBM,  an  actual  interceptor  of  some  kind  is  required. 
However,  since  the  objective  of  this  model  is  to  determine  the  extent  of  damage  to  a  target 
region  caused  by  TBM  fragments,  given  that  an  intercept  has  occurred,  modeling  of  the 
interceptor  is  omitted.  Therefore,  to  simulate  the  effect  of  intercept,  the  TBM  breaks 
(spontaneously)  apart  into  a  random  number  of  fragments,  each  experiencing  a  minimal 
change  in  velocity  from  that  of  the  TBM  because  most  of  the  energy  created  by  the 
intercept  is  consumed  fragmenting  the  TBM.  Also,  the  accumulated  mass  of  all  the 
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fragments  generated  is  taken  to  be  close  to  that  of  the  TBM.  Lacking  any  data  to  validate, 
these  assumptions  are  as  good  as  possible.  Accepting  these  two  conditions,  minimal 
changes  in  velocity  and  mass,  it  is  plausible  that  the  intercept  represents  a  worst-case 
situation. 

Determining  the  exact  number  of  fragments,  and  the  size  of  each  fragment,  that 
result  from  intercepting  a  TBM  is  impossible  because  of  numerous  factors  including  angle 
of  impact  by  interceptor,  angle  of  descent  of  ballistic  missile,  impact  point  on  ballistic 
missile,  distance  between  interceptor  and  ballistic  missile  at  intercept,  and  strength  of 
TBM  materials.  Therefore,  estimates  of  the  random  number  of  fragments  and  of  the 
individual  fragment  characteristics,  i.e.,  mass,  cross-sectional  area,  and  initial  velocity,  are 
introduced  to  complete  the  model. 

The  actual  realized  numbers  of  fragments  and  the  individual  fragment  masses  that 
occur  in  one  simulation  replication  are  developed  using  a  two-step  operation.  In  the  first 
step,  a  random  number  is  drawn  from  the  Poisson  probability  distribution.  The  mean 
number  of  fragments  desired  for  this  model,  which  is  requested  from  the  user  at  the  start 
of  the  simulation,  is  used  in  generating  this  Poisson  random  number.  This  number  is  then 
divided  into  the  TBM's  mass  to  get  the  average  mass  of  a  fragment.  The  second  step 
involves  a  loop  which  generates  fragments  one  at  a  time.  This  loop  is  dependent  solely  on 
the  mass  of  the  fragments  and  continues  until  the  accumulated  mass  of  the  fragments  is 
greater  than  that  of  the  TBM.  Each  time  through  the  loop  a  fragment  mass  is  generated 
by  multiplying  the  average  fragment  mass  produced  in  the  first  step  by  an  Exponential 
random  variable  with  a  mean  of  one.  This  process  simulates  the  variability  between 
fragment  masses.  Each  fragment  mass  is  then  added  to  the  accumulated  total.  So  long  as 
the  total  is  less  than  the  TBM's  mass,  a  fragment  is  generated.  When  the  accumulated 
total  becomes  greater  than  the  TBM's  mass,  fragment  generation  ceases  and  the  number  of 
fragments  is  determined.  Of  course  modifications  of  the  above  procedure  can  be  made  in 
the  simulation  code:  both  Poisson  and  Exponential  distributions  are  provisional  illustrative 
choices. 
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The  remaining  characteristics,  cross-sectional  area  and  initial  velocity,  of  each 
fragment  are  also  calculated  during  each  run  through  the  loop.  Three  assumptions  were 
made  in  order  to  specify  the  cross-sectional  area  of  a  fragment.  First,  since  the  exact 
shape  of  each  fragment  generated  is  unknown,  this  model  proposes  that  "the" 

'  o'is-seciional  area  is  a  function  of  the  fragment's  mass.  Second,  it  assumes  that  each 
iragmeni,  at  the  instant  of  intercept,  takes  on  the  (approximate)  shape  of  a  sphere.  The 
last  assumption  is  that  the  mass  density  of  the  fragment  equals  that  of  the  TBM.  With 
these  assumptions,  each  fragment's  cross-sectional  area  is  calculated  using  simple 
geometric  equations: 


Volume  = 

Density 

(18) 

Area  =  7t  •  radius^ 

(19) 

Once  again,  because  of  the  lack  of  knowledge  concerning  the  way  in  which  fragments  are 
gp  ■  ''■ated  after  an  intercept  of  a  TBM,  a  presumption  as  to  how  the  model  should  handle 
fragment  velocities  was  necessary.  The  initial  velocity  of  each  fragment  is  fairly  easy  to 
calculate,  assuming  the  fragments  are  dispersed  in  a  Gaussian  manner,  i.e.,  the  TBM 
breaks  up  in  a  conical  fashion,  and  that  there  is  negligible  loss  of  velocity.  To  accomplish 
f'’.is  Gaussian  dispersion,  a  random  "kick"  was  applied  to  the  three  components  of  the 
1  BM's  velocity  at  intercept.  The  amount  of  this  "kick"  was  determined  by  a  Normal 
random  variable  with  a  mean  of  zero  and  standard  deviation  of  ten  meters/second.  This 
allows  each  fragment  to  have  a  different  initial  velocity  and,  when  combined  with  different 
masses  and  cross-sectional  areas,  a  different  flight  path. 

3.  Fragment  Flight  Paths 

The  fragment  flight  paths  are  calculated  in  the  same  manner  as  is  that  of  the  TBM. 
A  loop  is  initiated  at  intercept  for  each  fragment  with  its  initial  conditions  estimated  under 
the  above  assumption.  At  each  pass  through  the  loop  the  fragment's  velocity  is  calculated 
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using  equations  (12)  -  (14)  from  Appendix  A.  The  fragments  experience  a  significant 
reduction  in  their  velocities.  This  change  is  the  result  of  greater  air  resistance  being 
applied  to  the  fragments  mainly  because  of  the  reduction  in  their  respective  masses. 
Fragments  with  larger  masses  experience  smaller  reductions  in  velocity  than  do  those  with 
smaller  masses.  The  positions  of  the  fragments  along  their  paths  are  then  updated  using 
equations  (1 5)  -  (17)  from  Appendix  A.  Prior  to  completing  each  pass  through  the  loop 
an  extra  step,  not  used  in  determining  the  flight  path  of  the  TBM,  is  implemented.  The 
extra  step  concerns  the  cross-sectional  area.  In  calculating  the  TBM's  flight  path,  which  in 
this  model  is  stable,  the  cross-sectional  area  remained  constant.  However,  it  is  believed 
that  the  flight  of  the  fragment  is  not  stable,  and  thus  the  cross-sectional  area  of  the 
fragment  does  not  remain  constant.  To  account  for  this  instability,  a  new  cross-sectional 
area  is  calculated  during  each  pass  of  the  loop.  This  is  accomplished  by  generating  a 
random  number  using  a  Normal  distribution  in  which  the  original  cross-sectional  area  acts 
as  the  mean,  and  the  standard  deviation  is  0.05  meters^.  One  might  speculate  that  the  drag 
coefficient,  C,  changes  along  with  the  changing  cross-sectional  area.  However,  that  case 
is  not  address  in  this  thesis.  One  thing  to  remember  is  that  this  procedure  is  simply  a 
means  to  simulate  the  unstable  flight,  and  does  not  represent  the  actual  way  in  which  a 
fragment  travels.  Upon  completion  of  a  fragment's  flight  path  the  loop  is  exited  and  the 
model  conducts  damage  assessment  by  locating  the  fragment's  hit  position  within  a  target 
subregion. 

4.  Target  Designation 

Defining  a  specific  target  in  the  model,  e  g.,  a  particular  populated  city  in  the 
middle  east,  is  too  restrictive  for  purposes  of  this  thesis.  The  objective  is  to  estimate  the 
damage  to  some  designated  target  region  caused  by  fragments  generated  from  an 
intercepted  TBM.  The  model  must  be  adaptable  for  different  targets.  To  accomplish  this, 
however,  one  restriction  is  necessary.  The  dimensions  of  the  target  region  must  be 
rectangular.  This  rectangular  region  is  then  divided  into  equal-sized  blocks,  or  subregions, 
to  allow  finer  detail  in  estimating  target  damage.  The  number  of  blocks  depends  on  the 
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desired  number  of  rows  and  columns  requested  from  the  user  at  the  start  of  the  simulation. 
The  generic  target  region  used  in  this  example  is  depicted  in  Figure  2.  The  exact 
center  of  the  rectangle  is  designated  as  the  aim  point.  This  rectangle  then  acts  as  a 
damage  template  which  is  placed  over  the  target  region  with  the  bottom  of  the  rectangle, 
blocks  13  through  16,  being  perpendicular  to  the  TBM  flight  path.  Figure  3  displays  a 
target  region  composed  of  blocks  6,  7,  10  and  11.  If  the  intended  target  is  not  rectangular 
in  shape  or  is  comprised  of  a  number  of  smaller  targets,  i.e.,  runways  at  a  coastal  airport 
or  piers  at  a  naval  port,  the  size  of  the  rectangle  is  increased  to  encompass  the  total  target 
area  and  only  the  blocks  containing  the  actual  targets  are  considered  areas  of  interest. 
Figure  4  depicts  piers  at  a  naval  port  composed  of  blocks  5,  6,  7,  9,  11,  13,  and  15.  The 
remaining  blocks  are  considered  areas  of  insignificance.  However,  a  target  that  is  located 
within  a  template  block  is  required  to  occupy  the  entire  block,  otherwise  the  damage 
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Figure  2.  Layout  of  generic  target  region  template 
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assessment  is  inaccurate  because  the  damage  suffered  is  not  uniform  throughout  a  block. 
This  encompassing  feature  of  the  target  rectangle  allows  the  target  to  be  of  any  size. 
However,  smaller  blocks  produce  more  accurate  damage  estimations  because  of  the  finer 
detail.  Assuming  that  all  targets  are  rectangle  is  fairly  unrealistic,  but  using  this 
assumption  for  damage  assessment  only,  allows  this  kind  of  model  to  be  used  for  just 
about  every  target  imaginable  that  fits  within  the  outlined  target  region. 

5.  Damage  Assessment 

This  model  assesses  damage  by  three  measures:  the  mean  number  of  fragments  to 
impact  the  target  subregions,  the  mean  total  mass  to  impact  the  target  subregions,  and  the 
mean  total  kinetic  energy  to  impact  the  target  subregions.  The  final  positions  of  each 
fragment  generated  from  a  single  intercept  are  compared  against  the  location  of  the 
designated  target  region  through  a  series  of  if-then  statements  to  determine  the  extent  of 
damage.  A  realistic  way  of  determining  damage  with  this  damage  template  is  to  assume 
that  the  entire  block  is  a  target,  i.e.,  there  are  no  insignificant  portions  of  the  target  block. 
Since  this  model  deals  with  a  generic  target  region,  a  generalization  is  made  about  the 
damage  that  occurs  to  a  target  block.  If  a  target  occupies  an  entire  block  then  the  damage 
that  occurs  to  that  block  is  treated  as  one  amount  even  though  fragments  may  impact  at 
several  different  points  within  the  block.  If  it  is  determined  that  the  fragment  impacted  the 
target  region,  the  estimated  damage  totals  of  the  target  block  which  contains  that 
fragment's  final  position  are  updated.  The  counter  for  the  number  of  fragments  to  impact 
is  incremented,  the  accumulator  for  the  total  mass  damage  has  the  mass  of  the  fragment 
added  to  it,  and  the  accumulator  for  the  total  kinetic  energy  damage  has  the  kinetic  energy 
of  the  fragment  added  to  it.  The  kinetic  energy  associated  with  fragment  /  is  calculated  by 
the  following: 

KE,  =  (20) 
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where 


V  =  velocity  of  fragment  /  at  impact 
m,  =  mass  of  fragment  /  at  impact 


In  order  to  produce  mean  damage  values  with  small  standard  errors,  at  least  2000 
replications  of  the  model  are  necessary  At  the  completion  of  each  replication  the  totals 
for  each  block  are  recorded  to  be  used  in  calculating  the  means  and  the  standard  errors  of 
the  means  using  the  following  formulas; 


R 


D.b  = 

^J=\ 


(21) 


'  =  fS-  2(0, i(/)  - 

ti-  1  y=l 


(22) 


where 


(23) 


Dib(j)  =  type  i  damage  value  for  block  b  in  the  jth  simulated  replication 
Dib  =  mean  type  i  damage  value  for  block  b 
R  =  number  of  replications 
s  =  sample  variance  of  damage  value 
SE  =  standard  error  of  the  mean 


In  addition  to  estimating  the  extent  of  damage  to  the  target  region,  the  model  also 
calculates  the  amount  of  damage  that  occurs  to  the  target's  surrounding  areas.  This  thesis 
considers  these  areas,  "miss"  areas.  However,  analysis  is  conducted  in  order  to  know  the 
location  of  the  remaining  fragment  damage  in  case  less  significant  targets  are  located  in 
these  areas.  Figure  5  displays  these  areas.  Since  the  launch  point  could  be  from  any 
direction,  assigning  the  surrounding  areas  with  a  direction,  i.e.,  north  of  the  target,  is  not 
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implemented.  Instead,  the  model  determines  if  the  damage  is  past  or  over  the  target,  or  is 
short  of  it,  or  right  or  left  of  it. 


Past  &  Left 

Past 

Past  &  Right 

Left 

TARGET 

Right 

Short  &  Left 

Short 

Short  &  Right 

Figure  5.  Plot  of  surrounding  areas  to  the  target  region. 


6.  Structure  of  a  Test 

Many  assumptions  were  made  in  this  model  concerning  the  fragmentation  of  a 
TBM  upon  interception.  The  reason  for  this  is  the  lack  of  available  data  on  both  the  true 
fragmentation  and  the  characteristics  of  the  resulting  fragments.  One  possible  way  of 
validating  this  model  is  to  structure  a  test  that  replicates  an  intercepted  TBM  and  the 
resulting  fragments.  After  completion  of  such  a  test,  a  distribution  fitted  to  the  fragment 
impact  points  determines  if  the  assumptions  made  were  correct  leading  to  a  valid  model. 
An  example  of  such  a  test,  on  a  small  and  inexpensive  scale,  could  use  2.75"  rockets  with 
ball  bearings  and  a  radio-controlled  burster  warhead  as  test  vehicles.  At  predetermined 
points  along  the  rocket's  trajectory  the  radio-controlled  burster  is  detonated  releasing  the 
ball  bearings  to  travel  their  own  paths  towards  the  target  region,  or  in  this  case  the  witness 
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area.  This  area  is  covered  in  saran  wrap  enabling  exact  determination  of  the  ball  bearing 
impact  points.  The  positions  of  these  impact  points  are  then  fitted  to  a  probability 
distribution  to  see  if  the  model  assumptions  are  correct.  This  is  just  one  example  of  a 
possible  test  that  can  be  used  to  validate  the  proposed  model  in  this  thesis  since  exact  data 
is  not  available. 

B.  DESIGN  OF  A  SIMULATION 

A  series  of  runs  are  conducted  to  estimate  the  damage  that  occurs  to  subregions 
within  a  designated  target  region  from  fragments  generated  by  the  intercept  of  a  TBM  as 
the  two  launch  angles,  vertical  and  off-target,  are  varied  slightly  to  account  for  the 
different  impact  points  associated  with  the  TBM.  To  increase  the  coverage  of  this  model 
and  examine  other  aspects  of  the  damage  estimates,  three  different  intercept  ranges  are 
simulated.  Varying  the  intercept  range  allows  an  understanding  of  the  need  to  intercept 
the  TBM  as  soon  as  possible,  thereby  causing  less  damage  to  the  target  area.  The 
MODSIM  source  code  for  the  simulation  is  contained  in  Appendix  C. 
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IV.  ILLUSTRATIVE  EXAMPLE 


To  aid  in  the  comprehension  of  this  model  an  illustrative  example  is  provided.  The 
TBM  is  launched  towards  a  target  located  560  km  downrange,  The  initial  conditions  for 
the  TBM  described  earlier  are  used;  vertical  launch  angle  of  30.0  degrees,  off-target 
launch  angle  of  0.0  degrees,  launch  acceleration  of  0.023875  meters/seconds^  ,  cross- 
sectional  area  of  0.6082  meters" ,  boost  phase  time  interval  of  90  seconds,  and  missile 
mass  of  6000  kg  after  boost  phase.  These  parameters  achieve  a  direct  hit  of  the  aim  point 
if  the  TBM  is  not  intercepted.  The  model  output  is  highly  sensitive  to  changes  in  these 
numbers. 

The  additional  information  necessary  to  conduct  this  simulation  are  the  intercept 
range  measured  in  kilometers  downrange  from  the  aim  point,  an  expected  number  of 
fragments  generated,  and  the  length  and  width  of  the  target  measured  in  kilometers.  For 
this  example,  the  intercept  ranges  of  interest  from  the  aim  point  are  2,  3,  and  4  kilometers, 
the  expected  number  of  fragments  is  5  [Ref  8],  and  2  km  for  the  length  and  width  of  the 
target,  corresponding  to  a  square  target  area  of  4  km^.  The  target  region  is  divided  into 
16  equal  blocks,  or  square  subregions,  each  of  area  0.25  km" .  Table  1  lists  the  estimates 
of  the  expected  number  of  fragments  to  impact  each  block,  the  expected  total  mass  to 
impact  each  block  and  the  expected  kinetic  energy  to  impact  each  block,  along  with  the 
standard  errors  associated  with  each  of  the  estimates  for  an  intercept  of  2  km.  Plots  of 
these  results  are  displayed  in  Figures  6  through  1 1 .  Table  2  show  the  results  for  an 
intercept  of  3  km,  and  the  plots  of  these  results  are  shown  in  Figures  12  through  17. 

Table  3  shows  the  results  for  an  intercept  of  4  km,  and  the  plots  of  these  results  are 
displayed  in  Figures  18  through  23.  The  remaining  vertical  launch  angles  examined  in  this 
thesis  are  the  set  (29.8° ,  29.9°,  30.1°,  30.2° ,  30.3° ).  The  remaining  ofF-target  launch 
angles  are  the  set  (  0.025°,  0.05°,  0.075°,  0.10°,  0.15°).  The  oflf-target  launch  angles 
represent  either  side  of  the  target  line.  The  model  was  run  for  the  right  side  of  the  target 
line.  To  get  the  data  for  angles  to  the  left  just  transpose  the  data  to  the  corresponding 
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blocks.  All  possible  combinations  of  these  angles  were  analyzed  at  the  three  intercepts 
and  the  results  are  contained  in  Appendix  B. 


Table  1 .  List  of  estimates  of  expected  values  and  corresponding  standard  errors  of 
the  estimates  for  number  of  fragment  impacts,  total  mass  impact  and  total 
kinetic  energy  impact  to  the  target  area  resulting  from  a  2  km  intercept. 
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Figure  6.  Mean  number  of  fragments  to  impact  the  target  which  is  divided  into  16 
equal  blocks  and  its  surrounding  areas.  This  data  is  generated  using 
launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  off-target  and 
an  intercept  range  of  2  km. 


X=  Downrange  Distance 
Y=  Off-Target  Distance 


Figure  7.  3-D  plot  of  mean  number  of  fragment  impacts  per  target  block  utilizing 
launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  off-target  and 
an  intercept  range  of  2  km. 
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.  ..ire  8.  Mean  total  mass  (kg)  to  impact  the  target  divided  into  16  equal  blocks  and 
its  surrounding  areas.  This  data  is  generated  using  launch  angles  of  30.0 
degrees  vertical  and  0.0  degrees  off-target  and  an  intercept  range  of  2  km. 


X=  Downrange  Distance 
Y=  Off-Target  Distance 


1258.15  kg 


Figure  9.  3-D  plot  of  mean  total  mass  (kg)  impact  per  target  block  utilizing  launch 
angles  of  30.0  degrees  vertical  and  0.0  degrees  off-target  and  an  intercept 
range  of  2  km. 


28 


Passed  & 

Passed 

Passed  & 

Left 

Right 

0 

0 

0 

0 

0 

0 

0 

Left 

0 

9.036 

9.036 

0 

Right 

0 

0 

18.81 

18.81 

0 

0 

0 

1.624 

1.624 

0 

Short  & 

Short 

Short  & 

Left 

Right 

0 

546.039 

0 

Figure  10.  Mean  total  kinetic  energy  (J)  to  impact  the  target  divided  into  16  equal 
blocks  and  its  surrounding  areas.  This  data  is  generated  using  launch 
angles  of  30,0  degrees  vertical  and  0.0  degrees  off-target  and  an 
intercept  range  of  2  km. 
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Figure  11.  3-D  plot  of  mean  kinetic  energy  (J)  impact  per  target  block  utilizing 

launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  off-target  and  an 
intercept  range  of  2  km. 


29 


Block  Number  Expected  Value  Standard  Error 

10-Frag 

0.273 

0.007 

10-Mass 

822.03  kg 

22.894  kg 

lO-KE 

12.182  J 

0.394  J 

1 1-Frag 

0.273 

0.007 

1 1-Mass 

822.03  kg 

22.894  kg 

n-KE 

12.182  J 

0.394  J 

0.471 

0.011 

14-Mass 

712.543  kg 

16.144  kg 

14-KE 

5.749  J 

0.134  J 

15-Frag 

0.471 

0.011 

15 -Mass 

712.543  kg 

16.144  kg 

15-KE 

5.749  J 

0.134  J 

Short-Frag 

3.137 

0.061 

29.438  kg 

Short-K  E 

510.411  J 

11.848  J 

Table  2.  List  of  estimates  of  expected  values  and  corresponding  standard  errors  of 
the  estimates  for  number  of  fragment  impacts,  total  mass  impact  and  total 
kinetic  energy  impact  to  the  target  area  resulting  from  a  3  km  intercept. 
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Figure  12.  Mean  number  of  fragments  to  impact  the  target  which  is  divided  into  16 
equal  blocks  and  its  surrounding  areas.  This  data  is  generated  using 
launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  off-target  and 
an  intercept  range  of  3  km. 


X=  Downrange  Distance 


Figure  13.  3-D  plot  of  mean  number  of  fragment  impacts  per  target  block  utilizing 
launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  off-target  and 
an  intercept  range  of  3  km. 
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Figure  14.  Mean  total  mass  (kg)  to  impact  the  target  divided  into  16  equal  blocks  and 
its  surrounding  areas.  This  data  is  generated  using  launch  angles  of  30.0 
degrees  vertical  and  0.0  degrees  ofF-target  and  an  intercept  range  of  3  km. 


X=  Downrange  Distance 


Y=  Off-Target  Distance 


Z=  Mean  Mass  Impact 


Figure  15.  3-D  plot  of  mean  total  mass  (kg)  impact  per  target  block  utilizing  launch 
angles  of  30,0  degrees  vertical  and  0.0  degrees  ofF-target  and  an  intercept 
range  of  3  km. 
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Figure  16.  Mean  total  kinetic  energy  (J)  to  impact  the  target  divided  into  16  equal 
blocks  and  its  surrounding  areas.  This  data  is  generated  using  launch 
angles  of  30.0  degrees  vertical  and  0.0  degrees  ofF-target  and  an 
intercept  range  of  3  km. 


X=  Downrange  Distance 


Y=  Off-Target  Distance 
Z=  Mean  KE  Impact 
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Figure  17.  3-D  plot  of  mean  kinetic  energy  (J)  impact  per  target  block  utilizing 

launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  ofif-target  and  an 
intercept  range  of  3  km. 
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Block  Number  Expected  Value  Standard  Error 

lO-Frag 

0.017 

0.002 

1  O’Mass 

88.549  kg 

10.683  kg 

lO’KE 

1.435  J 

0.174  J 

1  i-Frag 

0.017 

0.002 

1 1  “Mass 

88,549  kg 

10.683  kg 

11-KE 

1.435  J 

0.174  J 

14  “Frag 

0.259 

0.007 

14-Mass 

740.46  kg 

21.234  kg 

H-KE 

7.956  J 

0.241  J 

15 -Frag 

0.259 

0.007 

15 -Mass 

740.46  kg 

21.234  kg 

15-KE 

7.956  J 

0.241  J 

Shon-Frag 

4.073 

0.065 

Short-Mass 

2928.78  kg 

42.802  kg 

Short-K  E 

460.854  J 

10.895  J 

Table  3.  List  of  estimates  of  expected  values  and  corresponding  standard  errors  of 
the  estimates  for  number  of  fragment  impacts,  total  mass  impact  and  total 
kinetic  energy  impact  to  the  target  area  resulting  from  a  4  km  intercept. 
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Figure  1 8.  Mean  number  of  fragments  to  impact  the  target  which  is  divided  into 
1 6  equal  blocks  and  its  surrounding  areas.  This  data  is  generated  using 
launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  ofF-target  and 
an  intercept  range  of  4  km. 


X=  Downrange  Distance 


Y=  Off-Target  Distance 
Z=  Mean  Number  of  Hits 
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Figure  19.  3-D  plot  of  mean  number  of  fragment  impacts  per  target  block  utilizing 
launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  ofF-target  and 
an  intercept  range  of  4  km. 
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Figure  20.  Mean  total  mass  (kg)  to  impact  the  target  divided  into  16  equal  blocks 


and  its  surrounding  areas.  This  data  is  generated  using  launch  angles  of 
30.0  degrees  vertical  and  0.0  degrees  off-target  and  an  intercept  range 
of  4  km. 


X=  Downrange  Distance 


Y=  Off-Target  Distance 


Z=  Mean  Mass  Impact 
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Figure  21.  3-D  plot  of  mean  total  mass  (kg)  impact  per  target  block  utilizing 

launch  angles  of  30.0  degrees  vertical  and  0.0  degrees  off-target  and 
an  intercept  range  of  4  km. 
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Figure  22.  Mean  total  kinetic  energy  (J)  to  impact  the  target  divided  into  16  equal 
blocks  and  its  surrounding  areas.  This  data  is  generated  using  launch 
angles  of  30.0  degrees  vertical  and  0.0  degrees  ofF-target  and  an 
intercept  range  of  4  km. 


X=  Downrange  Distance 
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Figure  23.  3-D  plot  of  mean  kinetic  energy  (J)  impact  per  target  block  utilizing  launch 
angles  of  30.0  degrees  vertical  and  0.0  degrees  ofF-target  and  an 
intercept  range  of  4  km. 
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To  aid  in  data  interpretation,  it  is  assumed  that  the  target  damage  template 
encompasses  the  entire  target  region,  meaning  each  block  is  a  block  of  interest.  Also, 
damage  calculations  were  conducted  for  the  surrounding  areas  in  order  to  distinguish 
where  the  remaining  fragments  fall.  An  interesting  result  is  that  the  total  number  of 
fragment  impacts,  including  the  surrounding  areas,  experienced  in  the  three  intercept 
ranges,  as  shown  in  Figures  4,  6,  and  10,  are  equal.  This  is  expected  since  2000 
replications  are  run  and  all  impacts  are  tabulated.  The  total  masses  are  also  the  same  for 
the  three  intercept  ranges.  However,  as  expected,  the  total  kinetic  energies  decrease  as 
intercept  range  increases  because  the  terminal  velocities  of  the  fragments  are  smaller. 

Also  of  note  is  the  symmetry  experienced  around  the  target  region's  centerline  when  the 
otT-target  launch  angle  of  0.0”  is  used.  This  result  is  expected  when  taking  a  large  enough 
sample.  However,  with  only  2000  replications  a  small  degree  of  sampling  error  remained. 
To  account  for  this  error,  the  damage  values  that  resulted  for  the  blocks  on  either  side  of 
the  target  line  were  averaged  together  and  this  result  was  applied  to  both  blocks.  This 
procedure  is  justified  because  of  the  symmetry;  it  would  not  be  necessary  if  larger  samples 
were  taken  but  without  doing  so  the  computing  time  required  would  be  increased 
unnecessarily.  Any  other  off-target  launch  angle  does  not  produce  symmetry  around  the 
centerline  as  shown  in  Figure  24.  The  total  number  of  fragment  impacts  is  still  the  same 
with  this  new  off-target  launch  angle,  but  the  fragment  impact  positions  have  shifted  over 
to  the  right,  as  anticipated.  Increasing  this  angle  more  causes  a  further  shift  to  the  right  of 
the  fragment  impact  positions  up  to  a  point  where  no  fragments  impact  the  target  region. 

Determining  the  actual  resulting  damage  if  the  TBM  was  left  alone  to  impact  the 
target  region  at  the  aim  point  is  very  complicated.  Knowledge  of  any  structures  contained 
within  the  region  would  be  necessary,  along  with  the  manner  in  which  the  blast  from  the 
warhead  affects  these  structures.  To  simplify  this  problem,  the  same  damage  estimates 
used  for  the  fragments  are  employed.  Because  four  blocks  surround  the  aim  point  when 
launch  angles  of  30.0°  vertical  and  0.0°  off-target  are  used,  the  resulting  damage  is  spread 
equally  among  the  four  blocks  as  shown  in  Figures  25  and  26. 
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Figure  24.  Mean  number  of  fragments  to  impact  the  target  which  is  divided  into  16 
equal  blocks  and  its  surrounding  areas.  This  data  is  generated  using 
launch  angles  of  30.0  degrees  vertical  and  0.025  degrees  off-target  and 
an  intercept  range  of  2  km. 
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Figure  25.  Mean  total  mass  to  impact  target  region  when  TBM  is  not 
intercepted  and  aim  point  is  impacted. 


39 


Figure  26.  Mean  total  kinetic  energy  to  impact  target  region  when  TBM  is  not 
intercepted  and  aim  point  is  impacted 
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V.  DISCUSSION  AND  CONCLUSIONS 


The  estimates  of  damage  occurring  to  the  target  region  vary  depending  on  three 
factors  The  first  is  the  vertical  launch  angle.  As  this  angle  is  increased  the  number  of 
fragments  that  impact  the  target  region  increases,  along  with  the  respective  masses,  since 
the  TBM  is  at  a  higher  altitude  at  intercept  than  when  a  smaller  angle  is  employed.  This 
is  shown  in  Figures  6  and  27.  The  higher  altitude  produces  a  narrower  dispersion  pattern 
which  concentrates  the  fragments  impacts  closer  to  the  target.  However,  the  kinetic 
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Figure  27.  Mean  number  of  fragments  to  impact  the  target  which  is  divided  into  16 
equal  blocks  and  its  surrounding  areas.  This  data  is  generated  using 
launch  angles  of  30. 1  degrees  vertical  and  0.0  degrees  off-target  and 
an  intercept  range  of  2  km. 

energies  are  smaller  because  the  fragments  fly  longer  than  when  intercepted  at  a  lower 
altitude  thus  allowing  air  resistance  more  time  to  reduce  the  velocities.  The  second  factor 
is  the  off-target  launch  angle.  As  this  angle  is  increased,  i.e.,  moved  to  the  right  of  target 
line,  the  extent  of  damage  to  the  target  region  shifts  to  the  right  up  to  a  point  where  all  the 
fragments  miss  the  target  region.  A  shift  to  the  left  occurs  when  this  angle  is  decreased, 
i.e.,  moved  to  the  left  of  the  target  line.  Since  the  off-target  component  of  velocity  is 
treated  equally  in  both  cases,  moving  right  or  left  of  the  target  line,  the  extent  of  damage  is 


41 


equal,  however,  the  block  which  suffers  the  damage  is  transposed  to  the  opposite  side  of 
the  target  line.  As  shown  in  Appendix  B,  0. 15°  is  the  angle  which  shifts  everything  outside 
of  the  target  region.  The  shift  is  shown  in  Figures  6  and  24.  The  last  factor  which  causes 
variance  in  the  damage  estimates  is  the  intercept  range.  Intercepts  which  occur  close  to 
the  aim  point  generate  more  fragment  impacts  to  the  target  region  than  intercepts  which 
are  not  as  close.  The  difference  is  illustrated  in  Figures  6  and  18. 

These  three  factors,  taken  independently,  yield  fairly  distinct  distributions  of 
damage.  However,  when  examined  simultaneously,  these  factors  create  a  very  diverse 
damage  distribution.  Appendix  B  contains  the  layouts  of  all  possible  combinations  of 
these  three  factors  for  the  example  given  in  the  last  chapter.  The  only  exceptions  are  the 
off-target  angles  to  the  left  of  the  target  line.  The  extent  of  damage  is  the  same,  but  the 
blocks  in  which  the  damage  occurs  are  transposed  on  the  opposite  side  of  the  target  line. 
One  possible  approach  to  estimating  the  extent  of  damage  when  these  three  factors  are 
known  is  to  specify  a  response  function. 

Expected  Damage  =  /{Qot^Qv,  R)  (26) 


where 

©or  =  off-target  launch  angle 
0  F  =  vertical  launch  angle 

R  =  intercept  range  measured  in  km  downrange  from  the  target 

Such  a  function,  specified  on  physical  grounds,  could  be  fitted  to  data  obtained  at  various 
angles  and  ranges,  and  then  used  to  estimate  responses  at  others. 

The  conclusion  is  that  knowing  the  expected  damage  to  occur  to  a  target  region  if 
a  TBM  is  intercepted  is  helpful  in  making  a  decision  whether  to  intercept  it  or  not.  This  is 
obvious,  but  the  model  is  developed  to  illustrate  the  necessity  of  this  information.  This 
model  is  a  simple  analytical  attempt  to  produce  a  tentative  estimate.  At  this  time,  very 
little  is  known  concerning  the  dispersion  of  fragments  generated  by  an  intercept.  Because 
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of  this  lack  of  information,  a  crude  estimate  is  the  only  estimate  available.  When  the 
assumptions  made  in  this  thesis  are  determined  with  certainty,  a  more  detailed  version  of 
the  model,  replicating  the  exact  flight  paths  and  impact  points,  can  be  developed. 

Another  area  that  may  be  of  interest  is  determining  which  type  of  intercept 
approach  is  better,  a  blast-fragment  intercept  or  a  kinetic  energy  intercept.  A  blast- 
fragment  intercept  such  as  the  one  simulated  in  this  thesis  does  not  require  the  accuracy  of 
a  kinetic  energy  intercept  where  the  interceptor  must  actually  impact  the  TBM.  However, 
the  kinetic  energy  intercept  almost  completely  destroys  the  TBM  in  flight  resulting  in  very 
little  chance  of  collateral  damage.  The  technology  required  for  a  kinetic  energy 
interceptor  is  more  costly  than  that  of  the  blast-fragment  interceptor,  and  because  of  this 
high  cost  it  does  not  presently  exist.  The  Navy  is  experimenting  with  the  Light 
Exo-Atmospheric  Projectile  (LEAP),  a  kinetic  kill  vehicle  that  impacts  a  TBM  outside  of 
the  atmosphere.  If  this  experiment  is  successful  the  technology  may  be  applied  to  work 
within  the  atmosphere  on  intercepts  such  as  the  ones  proposed  in  this  thesis  eliminating 
the  need  for  an  estimate  of  fragment  damage. 
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APPENDIX  A. 


CALCULATION  OF  TBM  FLIGHT  PATH 

As  mentioned  in  the  text,  only  four  of  the  actual  characteristics  of  the  A1  Hussein 
are  used  in  this  model:  the  TBM's  mass,  the  TBM's  length,  the  TBM's  cross-sectional 
area,  and  the  TBM's  boost-phase  time  interval.  The  mass  of  the  TBM  at  launch  is  7000 
kg.  However,  at  the  end  of  boost  phase  the  mass  is  reduced  to  6000  kg  due  to  expended 
fuel.  The  method  for  compensating  for  expended  fuel  is  not  entirely  realistic  since  the 
mass  does  not  just  disappear  all  at  once.  However,  since  the  boost  phase  is  not  of  concern 
in  this  thesis,  and  since  the  assumption  has  been  made  that  straight-line  motion  during 
boost  was  unaffected  by  air  resistance  and  gravitation,  the  manner  in  which  missile  mass  is 
decreased  during  fuel  use  is  insignificant  so  long  as  the  mass  at  intercept  equals  the  mass 
at  the  completion  of  boost  phase.  The  length  of  the  TBM  is  12  meters.  The  diameter  of 
the  TBM  is  0.88  meters  at  its  widest  point  yielding  a  cross-sectional  area  of  0.6082 
meters'.  Finally,  determining  the  exact  boost-phase  time  interval  is  extremely  difficult  and 
requires  numerous  factors,  especially  the  fuel  burn  rate  of  the  TBM  which  depends  heavily 
upon  the  density  and  temperature  of  air.  Usually,  as  is  the  case  with  the  A1  Hussein,  the 
length  of  boost  phase  time  interval  is  given  by  an  approximation.  Because  of  this 
complexity  in  determining  the  exact  time  interval,  the  approximate  boost  phase  time 
interval  of  90  seconds  is  used  in  this  model  [Ref  2]. 

To  determine  the  remaining  parameters,  an  assumption  about  the  distance  to  the 
target  region  is  required;  it  is  that  the  center  of  the  target  region  designated  for  attack,  i.e., 
the  aim  point,  is  located  560  km  from  the  launch  point.  Having  knowledge  of  the  target 
distance  and  choosing  an  arbitrary  launch  angle  of  30  degrees  enables  the  launch 
acceleration  to  be  calculated  with  the  assumption  that  acceleration  was  constant 
throughout  the  entire  boost  phase.  The  calculated  launch  acceleration  is  0.023875 
meters/seconds^.  The  known  parameters,  together  with  these  parameters  calculated 
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under  the  above  assumptions,  are  the  only  data  required  to  model  the  launch  of  the  TBM 
and  begin  the  simulation.  Since  these  calculated  characteristics  are  not  the  precise 
characteristics  of  the  A1  Hussein  TBM,  the  flight  paths  modeled  in  this  simulation  do  not 
replicate  those  of  an  AJ  Hussein  exactly,  and  thus,  are  considered  generic  TBM  flight 
paths  usefully  close  to  the  actual  to  provide  insight.  All  of  the  above  numbers  can  be 
changed  in  the  software  at  the  will  of  an  analyst. 

TBM  Flight  Path  Calculations 

Prior  to  commencing  the  simulation,  the  user  is  prompted  for  three  items  that  aid 
in  the  construction  of  the  simulated  TBM  flight  path.  The  first  item  is  the  range  at  which 
intercept  is  going  to  occur,  measured  in  kilometers  downrange  from  the  aim  point.  The 
last  two  items,  the  vertical  and  the  off-target  launch  angles,  determine  the  location  of  the 
TBM's  impact  point  if  it  is  not  intercepted.  In  this  thesis  only  one  pair  of  angles,  30 
degrees  vertical  and  0  degrees  off-target,  result  in  a  direct  hit  of  the  aim  point.  All  other 
pairs  yield  different  impact  points  that  are  off-target-region  center.  Launching  a  TBM, 
even  one  with  state-of-the-art  technology,  at  a  target  560  km  away  results  in  some  level  of 
uncertainty  as  to  the  exact  impact  point.  Of  course,  better  technology  leads  to  smaller 
uncertainty.  Varying  the  vertical  and  off-target  launch  angles  slightly,  i.e.,  up  to  +/-  0.3 
degrees,  allows  parametric  sensitivity  analysis  to  be  conducted  in  this  model.  Such 
numbers  are  illustrative  and  can  be  altered  by  the  analyst. 

The  TBM  flight  is  broken  up  into  two  segments.  The  first  segment,  which  consists 
of  the  boost  phase,  is  fairly  easy  to  simulate  taking  into  account  the  assumption  of 
straight-line  motion,  during  which  air  resistance  and  gravity  have  negligible  effect.  Basic 
equations  of  physics  are  employed  to  determine  the  TBM's  position  and  velocity  at  the 
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completion  of  boost  phase.  The  linear  distance  traveled  by  the  TBM  during  boost  phase, 
D,  is  determined  from  a  side  calculation  where  gravity  and  thrust  are  taken  into  account. 
From  this  the  three-dimensional  position  of  the  TBM  at  the  completion  of  the  boost  phase 
is  calculated  using  the  following  equations: 

X  -  D  cos  {LA)  {OA)  (1) 

T  =  Z)  sin  {LA)  (2) 

Z  =  D  cos  {LA)  cos  {OA)  (3) 


where 


X  =  downrange  displacement  from  launch  point 
Y  =  altitude 

Z  =  off-target  displacement 
LA  =  vertical  launch  angle 
OA  =  off-target  angle 

The  TBM's  velocity  in  the  direction  tangent  to  flight  path  after  boost  phase,  Vq  ,  is  also 
determined  from  a  side  calculation  where  gravity  and  thrust  are  taken  into  account.  From 
this  the  components  of  velocity  are  calculated  in  the  following  manner 

Vx  =  Vo  cos  ( LA )  cos  ( OA )  (4) 


Vy  -  Vq  sin  ( LA ) 
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Vz  =  f  0  cos  ( LA )  sin  ( OA  )  (6) 

where 

Vv  =  downrange  component  of  velocity 
V^.  =  vertical  component  of  velocity 
=  off-target  component  of  velocity 

The  second  segment  consists  of  the  TBM  ballistic  flight  path.  The  purpose  for  its 
generation  is  to  enable  examination  of  various  intercept  points  so  as  to  determine  the 
lethality  of  the  fragments  created  by  those  intercepts  to  parts  of  the  target  region. 
Therefore,  it  is  necessary  to  be  able  to  calculate  the  positions  and  velocities  of  the  TBM 
throughout  its  entire  flight  path.  Initially,  an  attempt  is  made  to  replicate  the  exact  ballistic 
path  the  TBM  would  follow  through  the  atmosphere.  Knowing  that  gravity  and  air 
resistance  are  the  only  two  forces  which  affect  a  projectile's  flight,  an  expression  of  air 
resistance  is  sought  to  implement  in  a  closed-form  differential  equation  of  motion.  It  has 
been  determined  that  retardation  of  a  projectile  due  to  air  resistance  is  a  function  of  the 
projectile's  velocity  [Ref  3]. 

The  actual  flight  of  the  TBM  is  modeled  using  three-dimensional  equations  of 
motion  taking  gravity  and  atmospheric  density  into  consideration.  In  the  simulation,  a 
loop  is  created  which  continues  until  the  downrange  displacement  of  the  TBM  equals  the 
intercept  range  entered  by  the  user.  The  position  and  velocity  of  the  TBM,  calculated  by 
Equations  (1)  through  (6),  are  the  initial  conditions  prior  to  entering  the  loop.  For  ease  of 
computation,  an  interval  of  one  second  is  chosen  as  the  time  step  in  the  loop.  At  each 
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time  step  the  TBM's  terminal  velocity  at  time  t  is  calculated  by 


Vrit) 


I  2mg 

Cp(t)A 


where 


m  =  mass  of  the  TBM 

g  =  acceleration  due  to  gravity 

C  =  drag  coefficient  dependent  on  angle  of  attack 

p(t)  =  density  of  air  which  is  a  function  of  altitude  at  time  t 

A  =  cross-sectional  area  of  the  TBM 


Since  the  remaining  equations  are  all  functions  of  time,  the  reference  to  time,  i.e.,  (t),  is  a 
given  and  thus  removed  from  the  equations.  The  velocity  of  the  TBM  at  each  pass 
through  the  loop  is  updated  using 


y  =  (y}+  yj  +  r|)5 


From  these  two  calculations,  the  changes  in  each  component  of  velocity  are  calculated 
using  the  following  equations  of  motion 
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di^ 

dt 


g 


vv^ 


Y  g 


(10) 


til 


(11) 


From  these  equations  much  simpler  forms  for  finding  the  velocity  components  at  each 
( 1 -second)  time  step  are  defined  by 


-  - 2^ 


ri'y^CpA 

^  >  '  - r;;r- 


_  Cp^ 

2m 


(12) 

(13) 

(14) 


The  simple  one-step  differential  equation  solver  portrayed  above  can  be  replaced  by  far 
more  accurate  numerical  procedures,  such  as  4th-order  Runge-Kutta,  or  other 
methodologies  in  standard  package  programs  [Ref  4],  In  view  of  the  exploratory  nature 
of  this  investigation  such  complexity  has  not  been  introduced. 

Atmospheric  density,  p,  is  extremely  difficult  to  model.  The  Earth's  atmosphere  is 
primarily  composed  of  nitrogen  and  oxygen.  Solar  radiation  affects  the  dynamic 
properties  of  this  medium  by  constantly  changing  the  temperature,  pressure,  and  chemical 
constituents,  particulate  presence  and  electrical  properties  [Ref  5].  The  inability  to  model 
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atmospheric  density  creates  a  problem  when  trying  to  replicate  the  exact  flight  path  of  the 
TBM.  Therefore,  the  values  for  p  utilized  in  this  model  are  taken  from  the  U.S.  Standard 
•Atmosphere  Table  [Ref  6],  The  drag  coefficient,  C,  takes  into  account  the  relative 
contributions  of  viscous  and  form  resistances,  and  it  depends  on  the  nature  of  the  TBM 
(size,  shape,  and  irregularity  and  roughness  of  its  surface)  as  well  as  on  the  characteristics 
of  the  flow  of  air  over  the  TBM  [Ref  7],  Numerous  factors  are  involved  in  computing  C; 
a  value  of  1 .0  is  used  in  the  present  illustrative  simulation. 

The  position  of  the  TBM  is  also  required  at  every  step  in  the  loop.  Having 
calculated  the  velocity,  the  position  can  be  determined  rather  easily  by 


X, 

+ 

II 

(15) 

Vt 

II 

1 

t> 

+ 

(16) 

z, 

II 

•N 

7 

t> 

+ 

(17) 

The  loop  continues  until  the  downrange  displacement,  X, ,  equals  the  intercept  point 
determined  by  the  user.  When  the  intercept  point  is  reached  the  TBM  separates  into  a 
random  number  of  fragments. 
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Vertical  Launch  Angle  =  30.1°  Each  Block  Contains;  #  of  Fragments 

Off-Target  Launch  Angle  =  0.025°  Total  Mass 

Intercept  Range  =  4.0  km  Total  KE 
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P&L 

0 

Past 

0 

P&R 

0 

0 

0 

0 

0 

Left 

0 

0 

0.028 

140.454 

1.936 

0.012 

60.757 

0.834 

Right 

0 

0 

0 

0.241 

745.640 

7.670 

0.114 

356.967 

3.703 

0 

0 

0 

0.456 

843.933 

6.756 

0.168 

320.392 

2.600 

S  &  L 

0 

Short 

3.609 

2118.654 

634.174 

S&R 

0 

Vertical  Launch  Angle  =  30.  T 
Off-Target  Launch  Angle  =  0.05° 
Intercept  Range  =  4.0  km 

Each  Block  Contains:  #ofF 
Total 
Total 

P&L 

0 

Past 

0 

P&R 

0 

0 

0 

0 

0 

Left 

0 

0 

0 

0.040 

201.212 

2.770 

Right 

0 

0 

0 

0 

0.354 

1101.383 

11.362 

0 

0 

0 

0 

0.624 

1165.550 

9.366 

S&L 

0 

Short 

3.609 

2118.654 

634.174 

S&R 

0 

Vertical  Launch  Angle  =  30. 1° 
Off-Target  Launch  Angle  =  0.075° 
Intercept  Range  =  4.0  km 


Each  Block  Contains:  #  of  Fragments 
Total  Mass 
Total  KE 
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0 

0 

i 

0 

0 

0 

0 

0 

0 

3.54 
2044.002 
633.624 


0.543 

1011.775 

8.119 


Vertical  Launch  Angle  =  30.1° 
Off-Target  Launch  Angle  =  0.10° 
Intercept  Range  =  4.0  km 


Each  Block  Contains:  #  of  Fragments 
Total  Mass 
Total  KE 


Vertical  Launch  Angle  =  30.1° 
Off-Target  Launch  Angle  =  0.15° 
Intercept  Range  =  4.0  km 


1.014 

2463.733 

23.465 


S&R 

3.612 

2123.036 

634.207 


Each  Block  Contains:  #  of  Fragments 
Total  Mass 
Total  KE 
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?  Hi  L 

Past 

P&R 

0 

0 

0 

0 

0 

0 

0 

0.012 

0.012 

0 

66.104 

66.104 

0 

Left 

0.812 

0.812 

Right 

0 

0.098 

0.098 

0 

0 

355.413 

355.413 

0 

3.655 

3.655 

0.228 

0.228 

0 

529.662 

529.662 

0 

4.567 

4.567 

S&L 

Short 

S&R 

0 

3.950 

0 

2684.442 

471.931 

Vertical  Launch  Angle  =  30.2°  Each  Block  Contains:  #  of  Fragments 

Off-Target  Launch  Angle  =  0.0°  Total  Mass 


Intercept  Range  =  4.0  km 


Total  KE 


P&L 

0 

Past 

0 

P&R 

0 

0 

0 

0 

0 

Uft 

0 

0 

0.025 

132.207 

1.624 

0 

Right 

0 

0 

0 

0.196 

710.826 

7.310 

0 

0 

0 

0 

0.455 

1058.373 

9.126 

0 

S&L 

0 

Short 

3.951 

2685.392 

471.939 

S&R 

0 

Vertical  Launch  Angle  =  30.2° 
Off-Target  Launch  Angle  =  0.025° 
Intercept  Range  =  4.0  km 


Each  Block  Contains:  #  of  Fragments 
Total  Mass 
Total  KE 
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P&L 

Past 

P&R 

0 

0 

0 

0 

0 

0 

0 

0.018 

0.007 

0 

0 

94.449 

37.758 

Led 

1.161 

0.463 

Right 

0 

0.125 

0.071 

0 

0 

0 

454.326 

256.500 

4.672 

2.638 

! 

1 

i  0.303 

0.153 

0 

0 

:  703.010 

355.363 

i 

6.059 

3.067 

S&L 

Short 

S&R 

0 

3.951 

0 

2685.392 

471.939 

Vertical  Launch  Angle  =  30.2°  Each  Block  Contains:  #  of  Fragments 

Off-Target  Launch  Angle  =  0.05°  Total  Mass 

Intercept  Range  =  4.0  km  Total  KE 


Vertical  Launch  Angle  =  30.2°  Each  Block  Contains:  #  of  Fragments 
Off-Target  Launch  Angle  =  0.075°  Total  Mass 

Intercept  Range  =  4.0  km  Total  KE 
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P&L 

Past 

P&R 

0 

0 

0 

0 

0 

0 

0 

0.020 

0 

0 

0 

108.000 

Left 

1.328 

Right 

0 

0.147 

0.136 

0 

0 

0 

533.614 

5.484 

3.785 

0.372 

0 

0 

0 

864.897 

7.456 

S&L 

Short 

S&R 

0 

3.789 

0.163 

2477.783 

208.596 

470.465 

1.483 

Vertical  Launch  Angle  =  30.2°  Each  Block  Contains;  #  of  Fragments 

Off-Target  Launch  Angle  =  0. 10°  Total  Mass 


Intercept  Range  =  4.0  km 


Total  KE 


P&L 

0 

Past 

0 

P&R 

0 

0 

0 

0 

0 

Left 

0 

0 

0 

0 

Right 

0 

0 

0 

0 

0 

0.673 

1896.694 

18.023 

0 

0 

0 

0 

S&L 

0 

Short 

0 

S&R 

3.953 

2690.105 

471.977 

Vertical  Launch  Angle  =  30.2°  Each  Block  Contains:  #  of  Fragments 

Off-Target  Launch  Angle  =  0.15°  Total  Mass 

Intercept  Range  =  4.0  km  Total  KE 
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Vertical  Launch  Angle  =  30.3“  Each  Block  Contains:  #  of  Fragments 

Off-Target  Launch  Angle  =  0.0°  Total  Mass 


Intercept  Range  =  4.0  km 


Total  KE 


Intercept  Range  =  4.0  km 


Total  KE 
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P&L 

0 

Past 

0 

P&R 

0 

0 

0 

0.051 

239.215 

2.773 

0.027 

121.463 

1.385 

LciX 

0 

0 

0.184 

556.680 

5.307 

0.098 

298.826 

2.866 

Right 

0 

0 

0 

0.326 

657.157 

5.380 

0.157 

320.091 

2.634 

0 

0 

0 

0.396 

552.605 

3.989 

0.144 

198.627 

1.429 

S&L 

0 

Short 

3.246 

1642.134 

610.662 

S&R 

0 

Vertical  Launch  Angle  =  30.3°  Each  Block  Contains:  #  of  Fragments 

Off-Target  Launch  Angle  =  0.05°  Total  Mass 

Intercept  Range  =  4.0  km  Total  KE 


P&L 

0 

Past 

0 

P&R 

0 

0 

0 

0 

0.077 

360.679 

4.158 

Left 

0 

0 

0 

0.281 

855.505 

8.173 

Right 

0 

0 

0 

0 

0.482 

976.432 

8.007 

0 

0 

0 

0 

0.540 

752.049 

5.424 

S&L 

0 

Short 

3.246 

1642.134 

610.662 

S&R 

0 

Vertical  Launch  Angle  =  30.3° 
OfiF-Target  Launch  Angle  =  0.075° 
Intercept  Range  =  4.0  km 


Each  Block  Contains:  #  of  Fragments 
Total  Mass 
Total  KE 


105 


APPENDIX  B. 


RESULTS  OF  SIMULATION  WITH  VARIOUS  LAUNCH  ANGLES 


P  &  L 

Past 

P&R 

0 

0 

0 

0.060 

0 

0 

0 

282.068 

3.259 

0.217 

0 

0 

0 

656.479 

Lert 

6.254 

Right 

0 

0.393 

0.236 

0 

0 

0 

795.002 

548.573 

6.516 

4.954 

0.471 

0 

0 

0 

658.385 

4.753 

S&L 

Short 

S&R 

(J 

3.206 

0.044 

1602.535 

43.755 

610.074 

0.331 

Vertical  Launch  Angle  =  30.3°  Each  Block  Contains:  #  of  Fragments 

Off-Target  Launch  Angle  =  0.10°  Total  Mass 

Intercept  Range  =  4.0  km  Total  KE 


P&L 

0 

Past 

0 

P&R 

0 

0 

0 

0 

0 

Left 

0 

0 

0 

0 

Right 

0 

0 

0 

0 

0 

1.376 

2940.508 

25.736 

0 

0 

i 

0 

i 

0 

S&L 

0 

Short 

0 

S&R 

3.250 

1646.291 

610.405 

Vertical  Launch  Angle  -  30.3°  Each  Block  Contains:  #  of  Fragments 


Off-Target  Launch  Angle  =  0.15°  Total  Mass 

Intercept  Range  =  4.0  km  Total  KE 
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DEFINITION  MODULE  NewFrag; 

FROM  Ne\\TBM  IMPORT  G.  PositionT>T)e: 

CONST  C=  1.0; 

TYPE  TargetType  -  RECORD 
Rows;  INTEGER: 

Columns:  INTEGER: 

END  RECORD: 

IntArra>Type  =  ARRAY  INTEGER,  INTEGER  OF  INTEGER; 

RealArrayType  =  ARRAY  INTEGER,  INTEGER  OF  REAL; 

FragObj  =  OBJECT 
Mass  :  REAL; 

Area  :  REAL; 

Original  Area:  REAL; 

Velocin^  :  PositionType: 

Position  :  PositionType: 

Current Vel:  REAL: 

Damage  :  REAL: 

LengthofBlock:  REAL; 

Southwest:  PositionT>qDe; 

DX:  REAL; 

DZ:  REAL; 

TGT;  TargetType; 

ASK  METHOD  Objinit: 

ASK  METHOD  GetMass(IN  M:  REAL); 

ASK  METHOD  CalcArea: 

ASK  METHOD  CalcVelocity(IN  V:  PositionType); 

.^SK  METHOD  GetPosition(IN  P:  PositionType); 

ASK  METHOD  SetTargetRegion(IN  LE:  REAL;  IN  WI:  REAL;  IN  RO:  INTEGER;  IN  CO:  INTEGER;  IN  DI: 

REAL); 

ASK  METHOD  CalcDamage; 

ASK  METHOD  Travel; 

END  OBJECT; 

FragmentsType  =  ARRAY  INTEGER  OF  FragObj; 

VAR 

Fragment :  FragObj; 

Fragments  :  FragmentsType; 

FragCount :  IntArrayT\pe: 

MassTotal  :  RealArrayType; 

KETotal  :  RealArrayType; 

END  MODULE. 
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DEFINITION  MODULE  NewRunRcp: 
PROCEDURE  RiiiiRcplicaiions: 

END  MODULE. 


DEFINITION  MODULE  NewTBM: 

FROM  RandMod  IMPORT  RandomObj: 

FROM  NcwFrag  IMPORT  IiiiArrayType,  RealArrayType; 

CONST  G  -  .()()98: 

AcccI  =  ,023875: 

Area  =  6082: 

VVM  =  7000.0: 

MET  =  90.0: 
seedl  =  1: 
seed2  =  3 : 
seed3  =  5: 
seed4  =  7: 

Cl  =  1.0: 

TYPE  PositionTvpe  =  RECORD 
X;  REAL: 

Y;  REAL: 

Z;  REAL: 

END  RECORD: 

MissileObj  =  OBJECT 
Dist:  REAL: 

WarheadMass:  REAL: 

LA:  REAL: 

OA:  REAL: 

TBMPosition:  PositionT}pe; 

TBMVelocity^:  PositionType; 

Acceleration:  REAL; 

Pieces:  INTEGER: 

MolorBurnTime:  REAL: 

MissiieArea:  REAL: 

InierceptRange:  REAL: 

AfterMotorSpeed:  REAL; 

InterceptTime:  REAL; 

Altitude:  REAL; 

MeanFragments:  REAL; 

Length:  REAL: 

Width:  REAL: 

Rows:  INTEGER: 
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Columns:  INTEGER; 

TotalFragCount:  IntArrayType; 

TotalMD:  ReaiArrayType; 

TotalKE;  ReaiArrayType; 

TPLD:  REAL; 

TPLMD:  REAL; 

TPLKE:  REAL; 

TSLD:  REAL; 

TSLMD;  REAL; 

TSLKE:  REAL; 

TLD:  REAL; 

TLMD:  REAL; 

TLKE:  REAL; 

TPRD:  REAL; 

TPRMD:  REAL; 

TPRKE:  REAL; 

TSRD:  REAL; 

TSRMD;  REAL; 

TSRKE:  REAL; 

TRD:  REAL; 

TRMD:  REAL; 

TRKE:  REAL; 

TPD;  REAL; 

TPMD:  REAL; 

TPKE:  REAL; 

TSD:  REAL; 

TSMD:  REAL; 

TSKE:  REAL; 

ASK  METHOD  Objinit; 

ASK  METHOD  FlyTraj; 

ASK  METHOD  BreakUp; 

ASK  METHOD  GetInterceptRange(IN  R;  REAL); 

ASK  METHOD  GetMeanFragment(IN  MF;  REAL); 

ASK  METHOD  GetAngles(IN  VA;  REAL;  IN  OTA:  REAL); 

ASK  METHOD  GetAltitude(IN  ALT:  REAL); 

ASK  METHOD  GetTargetDimension(IN  L:  REAL;  IN  W;  REAL;  IN  R:  INTEGER;  IN  C:  INTEGER;  IN  D: 
REAL);  ASK  METHOD  ResetDamage(IN  Zero:  REAL); 

ASK  METHOD  GetDamage(IN  FC;  INTEGER;  IN  MT:  REAL;  IN  KET:  REAL;  IN  i:  INTEGER;  IN  j: 
INTEGER);  ASK  METHOD  GetPastLeftDamage(IN  DD:REAL;  IN  DMD:REAL;  IN  DKE:REAL); 

ASK  METHOD  GetShortLeftDamage(IN  DD:REAL;  IN  DMD:REAL;  IN  DKE:REAL); 

ASK  METHOD  GetLeftDamage(IN  DD:REAL;  IN  DMD:REAL;  IN  DKE:REAL); 

ASK  METHOD  GetPastRightDamage(IN  DD:REAL;  IN  DMD;REAL;  IN  DKE:REAL); 

ASK  METHOD  GetShortRightDamage(IN  DD:REAL;  IN  DMD:REAL;  IN  DKE;REAL); 

ASK  METHOD  GetRightDamage(IN  DD:REAL;  IN  DMD:REAL;  IN  DKE:REAL); 

ASK  METHOD  GetPastDamage(IN  DD:REAL;  IN  DMD:REAL;  IN  DKE;REAL); 

ASK  METHOD  GetShortDamage(IN  DD:REAL;  IN  DMD:REAL;  IN  DKE:REAL); 
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END  OBJECT. 

VAR  TBM:  MissilcObj; 
P.R.S.T:  RandomObj; 
END  MODULE. 


DEFINITION  MODULE  NewWriteLine; 
PROCEDURE  WriteLiiie(IN  Siring;  STRING); 
END  MODULE. 


IMPLEMENTATION  MODULE  NewFrag; 

FROM  RandMod  IMPORT  RandomObj; 

FROM  NewTBM  IMPORT  PositionType,  G,  TBM.R,S,T; 
FROM  MathMod  IMPORT  SQRT.POWER,  pi,CEIL; 
FROM  NewWriteLine  IMPORT  WriteLine; 


OBJECT  FragObj; 


ASK  METHOD  ObJInit; 
( 

\ . 

BEGIN 

NEW(Posilion); 

NEW(Velocity); 

NEW(TGT); 

NEW(SouthWest); 

END  METHOD; 


} 

} 


{ . } 

ASK  METHOD  GetMass(IN  M:  REAL); 

i . } 

BEGIN 
Mass  :=  M; 

END  METHOD; 


( . 

ASK  METHOD  CalcArea; 

{ . 

VAR 

Radius:  REAL; 


} 

} 
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BEGIN 

Radius  :=  ((3.0/(4.0*pi))  *  (0.001216425  *  Mass)); 
Area  pi  *  POWER(Radius,(2. 0/3.0)); 
OriginaiArea  :=  Area; 

END  METHOD: 

! . } 

ASK  METHOD  CalcVeIocity(IN  V:  PositionT>pe); 


kickx.kicia’.kickz  ;  REAL; 

BEGIN 

kickx  (ASK  T  TO  Normal(0.0,0.01)); 
kicky  :=  (ASK  T  TO  Normal(O.O.O.Ol)); 
kickz  (ASK  T  TO  Normal(0.0,0.01)); 

Velocity.X:=V.X  + kickx: 

Velocity. Y  V.Y  +  kick>- 
Velocity. Z  ;=  V.Z  +  kickz; 

END  METHOD: 

{ . } 

ASK  METHOD  GetPosition(IN  P:  PositionT>pe); 

{ . } 

BEGIN 

Position. X  :=  P.X; 

Position.Y  :=  P.Y; 

Positioii.Z  :=  P.Z: 

END  METHOD: 

{ . } 

ASK  METHOD  SetTargetRegion(IN  LE;  REAL;  IN  WI:  REAL;  IN  RO:  INTEGER;  IN  CO:  INTEGER;  IN  DI: 

REAL); 

{ . } 

BEGIN 

DX  :=  LE/FLOAT(RO); 

DZ  :=  WI/FLOAT(CO); 

South West.X  :=  DI  -  (0.5*FLOAT(RO)*DX); 

SouthWest.Z  ;=0.0  -  (0.5*FLOAT(CO)*DZ); 

TGT.Rows  :=  RO; 

TGT. Columns  :=  CO; 

NEW(FragCount.  1..RO,  L.CO); 

NEW(MassTotal.  1..RO,  L.CO); 

NEW(KETotaK  L.RO,  L.CO); 

END  METHOD; 
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I  I 

I . I 

ASK  METHOD  CalcDamage: 

{ . } 

VAR 

IJ.  PLDRDERD,SLDED,SRDXD.RD:  INTEGER; 

PLMDRLKERMDRKERRMDRRKE.SLMD,SLKE.SMD.SKE.  SRMD,SRKE,LMDXKE,RMD,RKE:  REAL; 
BEGIN 

I  :=  CElL({Posilion.X  -  SouthWest.X)/DX): 

J  CEIL((Position.Z  -  SouthWcst.Z)/DZ); 

IF  (I  >0)  AND(I<=  TGTRows) 

IF  (J  >  0)  AND  (J  <=  TGT.Columns) 

INC(FragCount[IJ]); 

MassTotai[lJ]  ;=  MassTotal[IJ]  +  Mass; 

KETotal[IJ]  KETotal[IJ]  +  (0.5*Mass*POWER(CurrentVeE2.0)); 

ASK  TBM  TO  GetDamage(FragCount[IJ],MassTotal[I,J],KETotal[IJ],I,J); 

END  IF; 

END  IF; 

IF  I  >  TGT.Rows 
IF  J  <=  0 
PLD  1; 

PLMD  Mass: 

PLKE  (0.5*Mass*POWER(CurrentVeI,2.0)); 

ASK  TBM  TO  GetPastLeftDaniage(FLOAT(PLD)XLMD,PLKE); 

END  IF; 

IF  (J  >  0)  AND  (J<=  TGT.Columns) 

PD  :=  1; 

PMD  :=  Mass; 

PKE  :=  (0.5*Mass*POWER(CurrentVeL2.0)); 

ASK  TBM  TO  GeLPastDamage(FLOAT(PD),PMDEKE); 

END  IF; 

IF  J  >  TGT.Columns 
PRD  :=  1: 

PRMD  ;=  Mass; 

PRKE  ;=  (0.5*Mass*POWER(CurrentVel.2.0)); 

ASK  TBM  TO  GetPastRightDamage(FLOAT(PRD), PRMD, PRKE); 

END  IF; 

END  IF; 

IF  I  <-  0 
IF  J  <=  0 
SLD  1; 

SLMD  ;=  Mass; 

SLKE  :=  (0.5*Mass*POWER(CurrentVel,2.0)); 

ASK  TBM  TO  GetShortLem:)amage(FLOAT(SLD),SLMD,SLKE); 

END  IF; 

IF  (J  >  0)  AND  (J  <=  TGT.Columns) 

SD  :=  1; 

SMD  Mass; 
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SKE  :=  (().5*Mass*POWER(CurrentVeI,2.0)); 

ASK  TBM  TO  GetShorLDamage(FLOAT(SD),SMD,SKE); 

END  IF: 

IF  J>TGT.Columns 
SRD  =  1: 

SRMD  :=  Mass: 

SRKE  :=  (0.5*Mass*POWER(CurrentVeK2.0)); 

ASK  TBM  TO  GetShortRightDamage(FLOAT(SRD), SRMD, SRKE); 
END  IF: 

END  IF: 

IF  (I  >  0)  AND  (I  <=  TGT.Rows) 

IF  J<=() 

LD  ;=  1: 

LMD  Mass: 

LKE  (0.5*Mass*POWER(CurrentVel,2.0)); 

ASK  TBM  TO  GetLeftDamage(FLOAT(LD).LMD,LKE); 

END  IF: 

IF  J  >=  TGT.Columns 
RD  I: 

RMD  :=  Mass: 

RKE  ;=  {0.5*Mass*POWER(CurrentVel,2.0)); 

ASK  TBM  TO  GetRightDamage(FLOAT(RD),RMD,RKE); 

END  IF: 

END  IF: 

END  METHOD: 

I . 

ASK  METHOD  Travel 


VAR 

L  rho.VTerm  :  REAL: 

BEGIN 
t  :=0.0: 

WHILE  (Position. Y  >  0.0)  AND  (Velocity.X  >  0.0) 
CASE  ROLrND(Position.Y) 

WHEN  0.. I: 
rho  :=  1. 1 1 17: 

WHEN  2: 
rho  1.0066: 

WHEN  3; 
rho  :=  0.90925: 

WHEN  4: 
rho  :=  0.81935; 

WHEN  5: 
rho  0.73643; 

WHEN  6: 
rho  :=  0.66011; 
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WHEN  7: 
rho  :=  0. 59002: 
WHEN  8: 
rho  :=  0.52579: 
WHEN  9: 
rho  :=  0.46706: 
WHEN  10: 
rho  :=  0  41351: 
WHEN  11; 
rho  :=  0.3648: 

WHEN  12: 
rho  :=  0.3 1 194: 
WHEN  13: 
rho  :=  0.2666: 

WHEN  14: 
rho  :=  0.22786: 
WHEN  15: 
rho  :=  0. 19475: 
WHEN  16: 
rho  :=  0.16647: 
WHEN  17: 
rho  ;=  0. 1423: 

WHEN  18: 
rho  :=  0. 12165; 
WHEN  19: 
rho  :=  0.104; 

WHEN  20: 
rho  :=  0.08991: 
WHEN  21. .25; 
rho  :=  0.04: 

WHEN  26.. 30: 
rho  :=  0.01841: 
WHEN  31.. 40; 
rho  :=  0.0039957: 
WHEN  41. .50: 
rho  :=  0.0010269: 
WHEN  51. .60: 
rho  :=  0.00030592; 
WHEN  61. .70: 
rho  0.000087535: 
WHEN  71. .80: 
rho  0.00001999: 
WHEN  81.. 150: 
rho  :=  0.0000004974: 
END  CASE: 
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CurrcnlVel  :=  SQRT(POWER(Velocity.X.2.0)  +  POWER(Velocity.Y,2.0)  +  POWER(Velocity.Z,2.0)); 
Vclociiy.X  :=  ({Vclocily.X*  1000.0)  -  ((CurrentVel’^1000,0*Velocin'.X*1000.0*rho*Area*C)/(2.0*Mass))) 
/lOOO.O; 

Velocity.Y  :=  ((Velocity. Y*1000.0)  -  (G*1000.0)  -  ((Current Vel*  1000.0* Velocity. Y*1000.0*rho*Area*C)/ 
(2.0*Mass)))/1000.0; 

Velocin^Z  :=  ((Velocity.Z*  1000.0)  -  ((CurrentVel*1000.0*Velocity.Z*1000.0*rho*Area*C)/(2.0*Mass))) 
/lOOO.O; 

Position.X  :=  Position.X  +  (Velocit}’.X); 

Posilion.Y  :=  Position. Y  +  (Velocity. Y); 

Position. Z  :=  Position. Z  +  (Velocity.Z); 

Area  :=  ASK  S  TO  Normal(OriginalArea.  0.05); 
l  ;=l+  1.0; 

END  WHILE; 

ASK  SELF  TO  CalcDamage; 

END  METHOD; 

END  OBJECT; 

END  MODULE. 


IMPLEMENTATION  MODULE  NewRunRep; 

FROM  NewTBM  IMPORT  TBM; 

FROM  SimMod  IMPORT  StartSimulation,ResetSimTime; 

FROM  NewWriteLine  IMPORT  WriteLine; 

FROM  MathMod  IMPORT  POWER,  SQRT; 

FROM  NewFrag  IMPORT  RealArrayType; 

PROCEDURE  RunReplications; 

VAR 

Z.p,q,c.d.ROW.COL:  INTEGER; 

APLD.APLMD,APLKE,ASLD,  ASLMD,ASLKE,ALD,ALMD,ALKE,APRD,APRMD,APRKE,ASRD:  REAL; 
ASRMD.ASRKE.ARD,ARMD,ARKE,APD,APMD,APKE,ASD,ASMD,ASKE,PLD2,PLMD2,PLKE2:REAL; 
SLD2.SLMD2  SLKE2,LD2,LMD2,LKE2,PRD2,PRMD2,PRKE2,SRD2,SRMD2,SRKE2,RD2;  REAL; 
RMD2.RKE2.PD2,PMD2,PKE2,SD2,SMD2,SKE2:  REAL; 
K,T,Range,VertAngle,OTAngle,MeanFragment,Len,Wid, Distance:  REAL; 

AvgF ragCount.  AvgMass,  AvgKE,  SQF ragCount,  SQMass,  SQKE:  RealArrayType; 

BEGIN 

NEW(TBM); 

UTPUTC'What  is  the  downrange  intercept  range  from  the  target?"); 

INPUT(Range); 

OUTPUTC'What  is  the  desired  vertical  launch  angle?"); 

INPUT(VertAngle); 

OUTPUTC'What  is  the  desired  off-target  launch  angle?"); 

INPUT(OTAngle); 

OUTPUTC'What  is  the  mean  number  of  fragments  desired?"); 

INPUT  (MeanFragment); 
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OUTPUT("What  is  the  length  of  the  target  in  KM?"); 

INPUT(Len); 

OUTPUTC'What  is  the  width  of  the  target  in  KM?"); 

INPUT(Wid); 

OUTPUTC'How  many  rows  of  blocks  in  the  target  region?"); 

INPUT(ROW); 

OUTPUTC'How  many  columns  of  blocks  in  the  target  region?"); 

INPUT(COL); 

OUTPUTC'What  is  the  distance  to  the  aimpoint  in  the  center  of  the  target?"); 

INPUT(Distance); 

ASK  TBM  TO  GetMeanFragment(MeanFragment); 

ASK  TBM  TO  GetAngles(VertAngle,OTAngle); 

ASK  TBM  TO  GetTargetDimension(Len.Wid.ROW.COL.Distance); 

ASK  TBM  TO  GetlntcrccptRange(Range); 

NEW(AvgFragCount.  1..ROW.  1..COL); 

NEW(AvgMass.  1..ROW,  1..COL); 

NEWCAvgKE.  1.  ROW.  1..COL); 

NEW(SQFragCoLint.  1..ROW.  L.COL); 

NEW(SQMass.  1..ROW.  1..COL); 

NEW(SQKE.  T.ROW.  1..COL); 

Z  :=  1; 

WHILE  Z<  101 
ASK  TBM  TO  FlyTraj; 

StartSimulation; 

FOR  c  :=  1  TO  ROW 
FOR  d  :=  1  TO  COL 

AvgFragCount[c.d]  :=  AvgFragCount[c,d]  +  FLOATfASK  TBM  TotalFragCount[c,d]); 
A\gMass(c.d]  ;=  AvgMass[c,d]  +  (ASK  TBM  TotalMD[c,d]); 

A\  gKEfc.d]  ;=  AvgKE[c.d]  +  (ASK  TBM  TotalKE[c,d]); 

SQFragCountIc.dJ  :=  SQFragCount[c,d]  +  POWER(FLOAT(ASK  TBM  TotaLFragCount[c,d]),2.0); 
SQMass[c,d]  :=  SQMass[c.d]  +  POWER(ASK  TBM  TotalMD[c.d],2.0); 

SQKE[c.d]  :=  SQKE[c.d]  +  POWER(ASK  TBM  TotalKE[c,d],2.0); 

END  FOR; 

END  FOR; 

APLD  :=  APLD  +  (ASK  TBM  TPLD); 

APLMD  :=  APLMD  +  (ASK  TBM  TPLMD); 

APLKE  :=  APLKE  +  (ASK  TBM  TPLKE); 

ASLD  :=  ASLD  +  (ASK  TBM  TSLD); 

ASLMD  :=  ASLMD  +  (ASK  TBM  TSLMD); 

ASLKE  :=  ASLKE  +  (ASK  TBM  TSLKE); 

ALD  :=  ALD  +  (ASK  TBM  TLD); 

ALMD  :=  ALMD  +  (ASK  TBM  TLMD); 

.ALKE  :=  ALKE  +  (ASK  TBM  TLKE); 

APRD  :=  APRD  +  (ASK  TBM  TPRD); 

APRMD  :=  APRMD  +  (ASK  TBM  TPRMD); 

APRKE  :=  APRKE  +  (ASK  TBM  TPRKE); 

ASRD  :=  ASRD  +  (ASK  TBM  TSRD); 

ASRMD  ;=  ASRMD  +  (ASK  TBM  TSRMD); 


116 


APPENDIX  C. 


MODSIM  SIMULATION  SOURCE  CODE 


ASRKE  :=  ASRKE  +  (ASK  TBM  TSRKE); 

ARD  :=  ARD  +  (ASK  TBM  TRD); 

ARMD  :=  ARMD  +  (ASK  TBM  TRMD); 

ARKE  :=  ARKE  +  (ASK  TBM  TRKE); 

APD  :=  APD  +  (ASK  TBM  TPD); 

APMD  :=  APMD  +  (ASK  TBM  TPMD); 

APKE  :=  APKE  +  (ASK  TBM  TPKE); 

ASD  :=  ASD  +  (ASK  TBM  TSD); 

ASMD  :=  ASMD  +  (ASK  TBM  TSMD); 

ASKE  =  ASKE  +  (ASK  TBM  TSKE); 

PLD2  :=  PLD2  +  POWER(ASK  TBM  TPLD.2.0); 

PLMD2  :=  PLMD2  +  POWER(ASK  TBM  TPLMD,2.0); 

PLKE2  :=  PLKE2  +  POWER(ASK  TBM  TPLKE.2.0); 

SLD2  :=  SLD2  +  POWER(ASK  TBM  TSLD.2.0); 

SLMD2  ;=  SLMD2  +  POWER(ASK  TBM  TSLMD.2.0); 

SLKE2  :=  SLKE2  +  POWER(ASK  TBM  TSLKE,2.0); 

LD2  ;=  LD2  +  POWER(ASK  TBM  TLD,2.0); 

LMD2  :=  LMD2  +  POWER(ASK  TBM  TLMD.2.0); 

LKE2  :=  LKE2  +  POWER(ASK  TBM  TLKE,2.0); 

PRD2  =  PRD2  +  POWER(ASK  TBM  TPRD,2.0); 

PRMD2  ■=  PRMD2  +  POWER(ASK  TBM  TPRMD,2.0); 

PRKE2  :=  PRKE2  +  POWER(ASK  TBM  TPRKE,2.0); 

SRD2  :=  SRD2  +  POWER(ASK  TBM  TSRD,2.0); 

SRMD2  :=  SRMD2  +  POWER(ASK  TBM  TSRMD,2.0); 

SRKE2  :=  SRKE2  +  POWER(ASK  TBM  TSRKE, 2.0); 

RD2  ;=  RD2  +  POWER(ASK  TBM  TRD,2.0); 

RMD2  :=  RMD2  +  POWER(ASK  TBM  TRMD,2.0); 

RKE2  :=  RKE2  +  POWER(ASK  TBM  TRKE,2.0); 

PD2  :=  PD2  +  POWER(ASK  TBM  TPD.2.0); 

PMD2  :=  PMD2  +  POWER(ASK  TBM  TPMD,2.0); 

PKE2  :=  PKE2  +  POWER(ASK  TBM  TPKE,2.0); 

SD2  :=  SD2  +  POWER(ASK  TBM  TSD,2.0); 

SMD2  :=  SMD2  +  POWER(ASK  TBM  TSMD,2.0); 

SKE2  :=  SKE2  +  POWER(ASK  TBM  TSKE, 2.0); 

ASK  TBM  TO  ResetDamage(O.O); 

ResetSimTime(O.O);  Z  ;=  Z  +  1; 

END  WHILE; 

K  :=  FLOAT(Z-l); 

T  :=K-  1.0; 

WriteLineC  "); 

WriteLineC’Intercept  Distance  is:  "  +  REALTOSTR(Range)  +  "  km"); 
WriieLineC'Vertical  Launch  Angle  is: "  +  REALTOSTR(VertAngle)  +  "  degrees"); 
WriteLineC’Off-Target  Launch  Angle  is;  "  +  REALTOSTR(OTAngle)  +  "  degrees"); 
WriieLineC'Intercept  Altitude  is:  "  +  REALTOSTR(ASK  TBM  Altitude)  +  "  km"); 
WriteLineC  "): 
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FOR  p  =  i  TO  ROW 
FOR  q  •=  1  TO  COL 

A\gFragCount|p.q]  :=  AvgFragCount|p.q]/K; 

AvgMasslp.q]  :=  AvgMass[p.q]/K; 

A\gKJE(p,ql  ;=  AvgKE|p.qJ/K; 

SQFragCoiintlp.qj  :=  SQRT(((SQFragCount[p.q]/T)  -  ((KCr)*POWER(AvgFragCount[p,q],2.0)))/K); 

SQMassIp.q]  :=  SQRT(((SQMass[p.ql/T)  -  ((m')*POWER(AvgMass[p,q],2.0)))/K); 

SQKEIp.q]  :=  SQRT(((SQKE[p.q]/T)  -  ((K/T)*POWER(AvgKE[p,q],2.0)))/K); 

WriicLincC’Block  "  -  INTTOSTR(p)  +  +  INTTOSTR(q)  +  "  MEAN  HITS:  "  + 

REALTOSTR(A\gFragCouiu|p.qJ)  +  "  Std  Error:  "  +  REALTOSTR(SQFragCount[p,q])); 

WriieLincC'Block  "  +  INTTOSTR(p)  +  +  INTTOSTR(q)  +  "  MEAN  MASS:  "  +  REALTOSTR(AvgMass[p,q]) 

-  "Std  Error:  "  +  REALTOSTR(SQMass(p,q])); 

WnicLincC'Block  "  -  INTTOSTR(p)  +  +  INTTOSTR(q)  +  "  MEAN  KE:  "  +  REALTOSTR(AvgKE[p,qJ) 

-  "  Std  Error:  "  +  REALTOSTR(SQKE[p.q])); 

END  FOR: 

END  FOR; 


APLD  :=  .APLD/K: 

APLMD  :=  APLMD/K; 

APLKE  :=  APLKE/K; 

ASLD  :=  ASLD/K; 

ASLMD  :=  ASLMD/K; 

ASLKE  :=  ASLKE/K; 

ALD  :=  ALD/K; 

ALMD  :=  ALMD/K; 

ALKE  ;=  ALKE/K; 

APRD  :=  APRD/K; 

APRMD  :=  APRMD/K; 

APRKE  ;=  APRKE/K: 

ASRD  :=  ASRD/K; 

ASRMD  :=  ASRMD/K; 

ASRKE  :=  ASRKE/K; 

ARD  :=  ARD/K: 

ARMD  ;=  ARMD/K: 

ARKE  :=  ARKE/K; 

APD  :=  APD/K; 

.APMD  :=  APMD/K; 

APKE  :=  APKE/K; 

ASD  :=  ASD/K; 

ASMD  :=  ASMD/K; 

ASKE  :=  ASKE/K; 

PLD2  ;=  SQRT((PLD2/T)-((K/T)*POWER(APLD,2.0)))/SQRT(K); 
PLMD2  :=  SQRT((PLMD2/T)-((K/r)*POWER(APLMD,2.0)))/SQRT(K); 
PLKE2  :=  SQRT((PLKE2/T)-((K/T)*POWER(APLKE,2.0)))/SQRT(K); 
SLD2  :=  SQRT((SLD2/T)-((KyT)*POWER(ASLD,2.0)))/SQRT(K); 
SLMD2  :=  SQRT((SLMD2/T)-((KTr)*POWER(ASLMD,2.0)))/SQRT(K); 
SLKE2  :=  SQRT((SLKE2/T)-((K/T)*POWER(ASLICE,2.0)))/SQRT(K); 
LD2  :=  SQRT((LD2Ar)-((KTr)*POWER(ALD,2.0)))/SQRT(K); 
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LMD2  :=  SQRT((LMD2/T)-((K/T)*POWER(ALMD,2.0)))/SQRT(K); 

LKE2  :=  SQRT((LKE2/T)-((KyT)*POWER(ALKE,2.0)))/SQRT(K); 

PRD2  :=  SQRT((PRD2/T)-((KyT)*POWER(APRD.2.0)))/SQRT(K); 

PRMD2  :=  SQRT((PRMD2/T)-((K/T)*POWER(APRMD.2.0)))/SQRT(K); 

PRKE2  :=  SQRT((PRKE2/T)-((KyT)*POWER(APRKE,2.0)))/SQRT(K); 

SRD2  :=  SQRT((SRD2/T)-((KyT)*POWER(ASRD,2.0)))/SQRT(K); 

SRMD2  •=  SQRT((SRMD2/T)-((KyT)*POWER(ASRMD,2.0)))/SQRT(K); 

SRKE2  :=  SQRT((SRKE2/T)-((K/T)*POWER(ASRKJE,2.0)))/SQRT(K); 

RD2  :=  SQRT((RD2/T)-((K/T)*POWER(ARD.2.0)))/SQRT(K); 

RMD2  :=  SQRT((RMD2/T)-((K/T)*POWER(ARMD,2,0)))/SQRT(K); 

RKE2  :=  SQRT((RKE2/T)-((K/T)*POWER(ARKE.2.0)))/SQRT(K); 

PD2  :=  SQRT((PD2/T)-((KyT)*POWER(APD,2.0)))/SQRT(K); 

PMD2  :=  SQRT((PMD2/T)-((K/T)*POWER(APMD.2.0)))/SQRT(K); 

PKE2  ;=  SQRT((PKE2/T)-((K/T)*POWER(APKE.2.0)))/SQRT(K); 

SD2  ;=  SQRT((SD2/T)-((K/T)*POWER(ASD.2.0)))/SQRT(K); 

SMD2  :=  SQRT((SMD2/T)-((K/T)*POWER(ASMD.2.0)))/SQRT(K); 

SKE2  :=  SQRT((SKE2/T)-((K/T)*POWER(ASKE.2.0)))/SQRT(K); 

WriieLineC'MISSES/LOCATIONS"); 

WriteLineC"’); 

WrileLineC'Mean  PastLeft  Misses:  "  +  REALTOSTR(APLD)  +  "  Std  Error:  "  +  REALTOSTR(PLD2)); 
WriteLineC’Mean  PastLeft  Mass  Damage:  "  +  REALTOSTR(APLMD)  +  "  Std  Error:  "  + 

REALTOSTR(PLMD2)); 

WriteLineC’Mean  PastLeft  KE  Damage:  "  +  REALTOSTR(APLKE)  +  "  Std  Error:  "  +  REALTOSTR(PLKE2)); 
WriteLineC'Mean  ShortLeft  Misses:  "  +  REALTOSTR(ASLD)  +  "  Std  Error:  "  +  REALTOSTR(SLD2)); 
WriteLineC'Mean  ShortLeft  Mass  Damage:"  +  REALTOSTR(ASLMD)  +  "  Std  Error:  "+REALT0STR(SLMD2)); 
WriteLineC'Mean  ShortLeft  KE  Damage:  "  +  REALTOSTR(ASLKE)  +  "  Std  Error:  "+  REALTOSTR(SLKE2)); 
WriteLineC’Mean  Left  Misses:  ’’  +  REALTOSTR(ALD)  +  "  Std  Error:  "  +  REALTOSTR(LD2)); 

WriteLineC’Mean  Left  Mass  Damage:  ’’  +  REALTOSTR(ALMD)  +  ’’  Std  Error:  ’’+  REALTOSTR(LMD2)); 
WriteLineC’Mean  Left  KE  Damage:  ’’  +  REALTOSTR(ALKE)  +  ’’  Std  Error:  ’’  +  REALTOSTR(LKE2)); 

WriteLineC’Mean  PastRight  Misses:  ’’  +  REALTOSTR(APRD)  +  ’’  Std  Error:  ’’  +  REALTOSTR(PRD2)); 
WriteLineC’Mean  PastRight  Mass  Damage:’’  +  REALTOSTR(APRMD)  +  "  Std  Error:  ’’+ 

REALTOSTR(PRMD2));  WriteLineC’Mean  PastRight  KE  Damage:  ’’  +  REALTOSTR(APRKE)  +  "  Std  Error:  ’’ 

+  REALTOSTR(PRKE2));  WriteLineC’Mean  ShortRight  Misses:  ’’  +  REALTOSTR(ASRD)  +  ’’  Std  Error:  ’’  + 
REALTOSTR(SRD2));  WriteLineC’Mean  ShortRight  Mass  Damage:’’  +  REALTOSTR(ASRMD)  +  ’’  Std  Error: 
’’+REALTOSTR(SRMD2));  WriteLineC’Mean  SkortRight  KE  Damage:  ’’  +  REALTOSTR(ASRKE)  +  ’’  Std  Error: 
’’  +  REALTOSTR(SRKE2));  WriteLineC’Mean  Right  Misses:  ’’  +  REALTOSTR(ARD)  +  ’’  Std  Error:  ’’  + 

REALTOSTR(RD2));  WriteLineC’Mean  Right  Mass  Damage:  ’’  +  REALTOSTR(ARMD)  +  ’’  Std  Error:  ’’  + 
REALTOSTR(RMD2));  WriteLineC’Mean  Right  BCE  Damage:  ’’  +  REALTOSTR(ARBCE)  +  ’’  Std  Error:  "  + 
REALTOSTR(RKE2));  WriteLineC’Mean  Past  Misses:  ’’  +  REALTOSTR(APD)  +  "  Std  Error:  ’’  + 

REALTOSTR(PD2));  WriteLineC’Mean  Past  Mass  Damage:  ’’  +  REALTOSTR(APMD)  +  "  Std  Error:  ’’  + 
REALTOSTR(PMD2));  WriteLineC’Mean  Past  BCE  Damage:  ’’  +  REALTOSTR(APBCE)  +  ’’  Std  Error:  ’’  + 

REALTOSTR(PBCE2));  WriteLineC’Mean  Short  Misses:  ’’  +  REALTOSTR(ASD)  +  ’’  Std  Error:  ’’  + 

REALTOSTR(SD2));  WriteLineC’Mean  Short  Mass  Damage:  ’’  +  REALTOSTR(ASMD)  +  ’’  Std  Error: "  + 
REALTOSTR(SMD2));  WriteLineC’Mean  Short  BCE  Damage:  ’’  +  REALTOSTR(ASBCE)  +  "  Std  Error:  "  + 
REALTOSTR(SBCE2)); 
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END  PROCEDURE: 
END  MODULE. 


IMPLEMENTATION  MODULE  NewTBM; 

FROM  NewWriteLine  IMPORT  WriteLine: 

FROM  SiniMod  IMPORT  SimTime: 

FROM  RaadMod  IMPORT  RandoniObj.FclchSeed: 

FROM  Ncw  Frag  IMPORT  FragObj.  FragmentsT\pe,  IntArrayType,  RealArravT\pe: 
FROM  MaihMod  IMPORT  COS,SIRSQRTPOWER,  pi; 

OBJECT  MissileObj: 

{ . } 

ASK  METHOD  Objinit: 

{ . } 

BEGIN 

NEW(TBM  Velocity); 

NBW(TBMPosition): 

NEW(P): 

NE\V(R): 

NEW(S): 

NEW(T): 

ASK  P  TO  SetSeed(FetchSeed(seedl)); 

ASK  R  TO  SetSeed(FetchSeed(seed2)); 

ASK  S  TO  SetSeed(FetchSeed(seed3)); 

ASK  T  TO  SetSeed(FetchSeed(seed4)); 

Acceleration  :=  Accel;  {kni/sec2} 

WarheadMass  :=  WM:  {kg} 

MissileArea  :=  Area;  {ni'^2} 

MotorBurnTime  :=  MBT  :  {secs} 

AfierMotorSpeed  :=  Acceleration  *  MotorBurnTime; 

END  METHOD: 

! . } 

ASK  METHOD  FlyTraj; 

{ . } 

VAR 

InterceptPt,  rho,  Current Vel:  REAL; 

BEGIN 

NEW(TotalFragCount.  L.Rows,  1.. Columns); 

NEW(TotalMD.  L.Rows.  1.. Columns); 

NEW(TotalK£.  L.Rows,  L. Columns); 

InterceptPt  :=  Dist  -  InterceptRange; 
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TBMPosition.X  :=  (Accel*0.5*MBT*MBT*COS(LA)*COS(OA)); 

TBMPosition.  Y  =  (Accel*0.5*MBT*MBT*SIN(LA)): 

TBMPosition.Z  :=  (Accel*().5*MBT*MBT*COS(LA)*SIN(OA)); 

TBMVclociiy.X  (Accei*MBT)*COS(LA)*COS(OA); 

TBMVclociiy.Y  ;=  (Accel*MBT)*SIN(LA); 

TBMVelocih  .Z  :=  (Accel*MBT)*COS(LA)*SIN(OA); 

WarheadMass  :=  6000.0; 

WHILE  FLOAT(ROUND(TBMPosition.X))  <  InterceptPt 

CurrentVel  :=  SQRT(POWER(TBMVelocity.X,2.0)+POWER(TBMVelocity.Y,2.0)  + 
POWER(TBMVelocity.Z.2.0)); 

IF  TBMPosition.  Y  >  0.0 
CASE  ROUND(TBMPosilion.Y) 

WHEN()..l; 
rho  =  1  1117. 

WHEN  2: 
rlio  1 .0066: 

WHEN  3; 
rho  0.90925: 

WHEN  4; 
rho  :=  0.81935; 

WHEN  5: 
rho  :=  0.73643; 

WHEN  6: 
rho  :=  0.66011: 

WHEN  7; 
rho  :=  0.59002: 

WHEN  8: 
rho  :=  0.52579: 

WHEN  9: 
rho  :=  0.46706: 

WHEN  10: 
rho  :=  0.41351; 

WHEN  11: 
rho  0.3648; 

WHEN  12: 
rho  :=  0.31194; 

WHEN  13: 
rho  :=  0.2666: 

WHEN  14: 
rho  :=  0.22786: 

WHEN  15: 
rho  :=  0.19475; 

WHEN  16: 
rho  :=  0.16647: 

WHEN  17: 
rho  :=  0.1423; 
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WHEN  18: 
rho  :=  0.12165: 

WHEN  19: 
rho  :=  0.104: 

WHEN  20: 
rho  :=  0.08991: 

WHEN  21. .25: 
rho  :=  0.04: 

WHEN  26.. 30: 
rho  :=  0.01841: 

WHEN  31. .40: 
rho  :=  0.0039957: 

WHEN4L.50: 
rho  :=  0.0010269: 

WHEN  51. .60: 
rho  :=  0.00030592: 

WHEN  61. .70: 
rho  :=  0.000087535: 

WHEN  71. .80: 
rho  0.00001999: 

WHEN  81. .200: 
rho  0.0000004974: 

END  CASE: 

END  IF; 

TBMVelocit>^X  :=  ((TBMVelocity.X*  1000.0)  -  ((Current Vel*  1000.0*  TBMVelocity.X*  1000.0* 
rho*MissileArea*Cl)  /(2.0*WarheadMass)))/1000.0; 

TBMVelocily.Y  :=  ((TBMVelocity.Y*  1000.0)  -  (0.0098*1000.0) « ((Current Vel*l 000.0 *TBM Velocity. Y* 
1000.0*rho*MissileArea*Cl)/(2.0*Warhea(iMass)))/1000.0: 

TBMVelociU-Z  :=  ((TBMVelociW-Z*  1000.0)  -  ((CurrentVel*1000.0*TBMVelocity.Z*  1000.0*  rho* 
MissileArea*  Cl)/(2.0*WarheadMass)))/1000.0: 

TBMPosition.X  :=  TBMPosition.X  +  TBMVelocity.X; 

TBMPosition.Y  :=  TBMPosition.Y  +  TBMVelocity.Y; 

TBMPosiiion.Z  :=  TBMPosition.Z  +  TBMVelocity.Z: 

END  WHILE: 

ASK  SELF  TO  BreakUp; 

ASK  SELF  TO  GetAltitude(TBMPosition.  Y); 

END  METHOD; 

{ . } 

ASK  METHOD  BreakUp; 

( . } 

VAR 

l. K:  INTEGER: 

m. TotalMass:  REAL: 

Fragment:  FragObj: 

Fragments:  FragmentsType; 
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BEGIN 

TotalMass  :=  0.0: 

K  =  ASK  R  TO  Poisson(ASK  SELF  MeanFragments); 

IF  K  <=  I 
K;-2: 

END  IF: 

NEW(Fragments,  L.5*K); 

I  :=  1; 

WHILE  TotalMass  <  WarheadMass 

m  :=  (WarheaclMass/FLOAT(K))*(ASK  P  TO  Exponential(LO)); 

TotalMass  :=  TotalMass  +  m: 

IF  m  >  70.0 

IF  TotalMass  <  WarheadMass 
NEW(FragmeiusllJ): 

ASK  Fragmcnts|I]  TO  GctMass(ni): 

ASK  Fragments|I]  TO  CalcArea: 

ASK  Fragments[I]  TO  CalcVelocity(TBMVelocity); 

ASK  Fragments[I]  TO  GetPosition(TBMPosition); 

ASK  Fragments[I]  TO  SetTargetRegion(Length, Width, Rows, Columns,  Dist); 
ASK  Fragments[I]  TO  Travel; 

I  :=I+  1: 

END  IF; 

END  IF: 

END  WHILE: 

END  METHOD: 


! . } 

ASK  METHOD  GetInterceptRange(IN  R:  REAL); 

{ . } 

BEGIN 

nterceptRange  :=  R; 

END  METHOD; 


{ . } 

ASK  METHOD  GetMeanFragment(IN  MF:  REAL); 


BEGIN 


} 


MeanFragments  :=  MF; 
END  METHOD; 


{ . } 

ASK  METHOD  GetAngles(IN  VA:  REAL;  IN  OTA:  REAL); 

{ . } 

BEGIN 

LA  :=  (VA*pi)/180.0; 

OA  :=  (OTA*pi)/180.0; 

END  METHOD; 
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APPENDIX  C. 


MODSIM  SIMULATION  SOURCE  CODE 


<  ) 

I . / 

ASK  METHOD  GetAltitude(IN  ALT:  REAL); 

{ . } 

BEGIN 

Allitudc  •=  ALT: 

END  METHOD: 

! . } 

ASK  METHOD  GciTargeiDimcnsion(IN  L:  REAL;  IN  W:  REAL;  IN  R:  INTEGER;  IN  C;  INTEGER;  IND 

REAL); 

»  > 

^ . f 

BEGIN 
Length  :=  L: 

Width  :=  W: 

Rows  :=  R: 

Columns  :=  C: 

Dist  :=  D: 

END  METHOD: 

f . } 

ASK  METHOD  GeiDamage(IN  FC;  INTEGER;  IN  MT:  REAL;  IN  KET;  REAL;  IN  i;  INTEGER;  IN  j; 

INTEGER); 

. . } 

BEGIN 

TotalFragCount[i,j]  :=  TotaIFragCount[i.j]  +  FC; 

TotalMD[i.j]  ;=  TotalMD[ij]  +  MT; 

TotalKEfiJ]  :=  TotalKE[ij]  +  KET; 

END  METHOD; 

i . } 

ASK  METHOD  GetPastLeftDamage(IN  DD:  REAL;  IN  DMD;  REAL;  IN  DKE;  REAL); 

{ . } 

BEGIN 

TPLD  ;=  TPLD  +  DD; 

TPLMD  :=  TPLMD  +  DMD; 

TPLKE  :=  TPLKE  +  DKE; 

END  METHOD; 

( . } 

ASK  METHOD  GetShortLeftDamage(IN  DD:  REAL;  IN  DMD:  REAL;  IN  DICE'  REAL)' 

{ . } 

BEGIN 

TSLD  :=  TSLD  +  DD; 

TSLMD  :=  TSLMD  +  DMD; 

TSLKE  :=  TSLKE  +  DKE; 

END  METHOD; 
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MODSIM  SIMULATION  SOURCE  CODE 


ASK  METHOD  GciLeflDamage(IN  DD 


BEGIN 

TLD  ;=  TLD  +  DD; 
TLMD  :=  TLMD  +  DMD; 
TLKE  :=  TLKE  +  DKE; 
END  METHOD; 


I . > 

ASK  METHOD  GelPaslRigluDamage(lN  DD:  REAL;  IN  DMD:  REAL;  IN  DKE;  REAL); 

•  ) 

I . f 

BEGIN 

TPRD  =  TPRD  +  DD; 

TPRMD  :=  TPRMD  +  DMD; 

TPRKE  :=  TPRKE  +  DICE; 

END  METHOD; 

{ . } 

ASK  METHOD  GetShortRightDamage(IN  DD;  REAL;  IN  DMD:  REAL;  IN  DKE;  REAL); 

{ . } 

BEGIN 

TSRD  :=  TSRD  +  DD; 

TSRMD  :=  TSRMD  +  DMD; 

TSRKE  :=  TSRKE  +  DKE; 

END  METHOD; 

{ . } 

ASK  METHOD  GetRightDamage(IN  DD:  REAL;  IN  DMD:  REAL;  IN  DKE;  REAL); 

{ . } 

BEGIN 

TRD  :=  TRD  +  DD; 

TRMD  :=  TRMD  +  DMD; 

TRKE  :=TRKE  +  DKE; 

END  METHOD; 


:  REAL:  IN  DMD:  REAL;  IN  DKE:  REAL); 


} 


ASK  METHOD  GetPastDamage(IN  DD:  REAL;  IN  DMD:  REAL;  IN  DKE:  REAL); 

{ . } 


BEGIN 


TPD  ;=  TPD  +  DD; 
TPMD  :=  TPMD  +  DMD; 
TPKE  :=  TPKE  +  DKE; 
END  METHOD; 
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MODSIM  SIMULATION  SOURCE  CODE 


{ . } 

ASK  METHOD  GclShortDaniagc(IN  DD:  REAL;  IN  DMD:  REAL;  IN  DKE:  REAL) 

‘  > 

< . r 

BEGIN 

TSD  ;=TSD  ^  DD: 

TSMD  ;=  TSMD  ^  DMD: 

TSNE  :=TSKE-t-DKE: 

END  METHOD: 


. i 

ASK  METHOD  ResetDamage(IN  Zero:  REAL); 

{ . } 

VAR 

k.n:  INTEGER: 

BEGIN 

FOR  k  I  TO  Rows 

FOR  n  I  TO  Columns 

TolaIFragCount[k,n]  :=  0: 

TotalMD[k,n]  Zero: 

TotalKE[k,n]  :=  Zero: 

END  FOR: 

END  FOR: 

TPLD  :=  Zero: 

TPLMD  :=  Zero: 

TPLKE  ;=  Zero: 

TSLD  :=  Zero: 

TSLMD  :=Zero: 

TSLKE  :=  Zero: 

TLD  :=  Zero: 

TLMD  :=  Zero: 

TLKE  :=  Zero: 

TPRD  :=  Zero: 

TPRMD  ;=  Zero: 

TPRKE  :=  Zero; 

TSRD  :=  Zero; 

TSRMD  :=  Zero; 

TSRKE  :=  Zero; 

TRD  :=  Zero; 

TRMD  :=  Zero: 

TRKE  ;=  Zero: 

TPD  :=  Zero; 

TPMD  :=  Zero: 

TPKE  :=  Zero; 

TSD  :=  Zero: 

TSMD  :=  Zero; 


APPENDIX  C. 


MODSIM  SIMULATION  SOURCE  CODE 


TSKE  :=  Zero: 
END  METHOD: 
END  OBJECT: 

END  MODULE. 


IMPLEMENTATION  MODULE  NewWriteLine: 

FROM  lOMod  IMPORT  FileUseType(Output): 
FROM  lOMod  IMPORT  SlreamObj: 

FROM  UliIMod  IMPORT  DateTime: 

FROM  SiniMod  IMPORT  SimTime; 

FROM  NewTBM  IMPORT  TBM: 

VAR 

DT  ;  STRING: 

TraceStream  :  StreamObj; 

PROCEDURE  WriteLine(IN  String:  STRING): 
BEGIN 

IF  (TraceStream  =  NILOBJ) 

NEW(TraceStream): 

ASK  TraceStream  TO  Open("sim _ .out”.  Output); 

DateTime(DT); 

ASK  TraceStream  TO  WriteString(DT); 

ASK  TraceStream  TO  WriteLn; 

ASK  TraceStream  TO  WriteLn;  E 
NDIF; 

ASK  TraceStream  TO  WriteReaI(SimTime(),  7,3); 
ASK  TraceStream  TO  WriteString("  —  "  +  String); 
ASK  TraceStream  TO  WriteLn: 

END  PROCEDURE: 

PROCEDURE  WriteLineCIose: 

BEGIN 

ASK  TraceStream  TO  Close; 

END  PROCEDURE; 

END  MODULE. 
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